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ABSTRACT 

The carbon-enhanced metal-poor (CEMP) stars constitute approximately one fifth of the metal-poor ([Fe/H] < —2) 
population but their origin is not well understood. The most widely accepted formation scenario, at least for the majority 
of CEMP stars which are also enriched in s-process elements, invokes mass-transfer of carbon-rich material from a 
thermally-pulsing asymptotic giant branch (TPAGB) primary star to a less massive main-sequence companion which is 
seen today. Recent studies explore the possibility that an initial mass function biased toward intermediate-mass stars is 
required to reproduce the observed CEMP fraction in stars with metallicity [Fe/H] < —2.5. These models also implicitly 
predict a large number of nitrogen-enhanced metal-poor (NEMP) stars which is not seen. In this paper we investigate 
whether the observed CEMP and NEMP to extremely metal-poor (EMP) ratios can be explained without invoking 
a change in the initial mass function. We construct binary-star populations in an attempt to reproduce the observed 
number and chemical abundance patterns of CEMP stars at a metallicity [Fe/H] ~ —2.3. Our binary-population models 
include synthetic nucleosynthesis in TPAGB stars and account for mass transfer and other forms of binary interaction. 
This approach allows us to explore uncertainties in the CEMP-star formation scenario by parameterization of uncertain 
input physics. In particular, we consider the uncertainty in the physics of third dredge up in the TPAGB primary, 
binary mass transfer and mixing in the secondary star. We confirm earlier findings that with current detailed TPAGB 
models, in which third dredge up is limited to stars more massive than about 1.25 M Q , the large observed CEMP 
fraction cannot be accounted for. We find that efficient third dredge up in low-mass (less than 1.25 Mq), low-metallicity 
stars may offer at least a partial explanation to the large observed CEMP fraction while remaining consistent with the 
small observed NEMP fraction. 

Key words. Stars: carbon - Stars: binaries - Stars: chemically peculiar - Galaxy: halo - Galaxy: stellar content - 
Nucleosynthesis 



1. Introduction 

One of the most interesting problems in modern stel- 
lar astronomy is to explain the existence of a popula- 
tion of carbon-enhanced metal-poor (CEMP 1 ) stars in 
the Galactic halo. The HK (Beers et ah, 1992) and 
Hamburg/ESO (Christlieb et al., 2001) surveys find a 
large number of CEMP stars among the metal-poor (EMP, 
[Fc/H] < —2) population, at a fraction around 20% (e.g. 
9 ± 2% Frebel et al., 2006; 21 ± 2% Lucatello et al., 2006; 
up to 30% from the SAGA database of Suda et al., 2008 - 
see below for details). 

The CEMP stars are subdivided into four groups de- 
pending on the presence or absence of the heavy elements 
barium and europium (see e.g. Beers & Christlieb, 2005). 
The most populous group consists of the s-process rich 
CEMP stars, the so-called CEMP-s stars (e.g. Aoki et al., 
2007) which display barium enhancements of [Ba/Fe] > 
+0.5. These account for about 80 per cent of all CEMP 
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1 [C/Fe] > 1, [Fe/H] < —2, where the logarithmic abundance 
ratio [X/Y] = log 10 (X/Y) - log lo (X / Y ). 



stars. There are also CEMP stars with r-process enhance- 
ments (the CEMP-r class) and some with both r- and 
s-process enhancements (CEMP-r+s, e.g. Jonsell et al., 
2006) . Finally, there is a class of CEMP stars which show no 
enhancement of neutron-capture elements. These are called 
the CEMP-no stars (Aoki et al., 2002). A detailed review 
of the various CEMP subgroups can be found in Masseron 
et al. (2009). 

A quantitative understanding of the origin of CEMP 
stars touches on many branches of stellar astronomy. The 
most likely formation mechanism for the s-process rich 
CEMP stars involves mass transfer in binary systems. 
Carbon-rich material from the TPAGB primary star pol- 
lutes the lower-mass main sequence secondary such that it 
becomes enriched in carbon and s-process elements. We ob- 
serve only the secondary today; the primary is an unseen 
white dwarf. Surveys of radial velocity shifts find that the 
binary fraction of CEMP-s stars is consistent with them 
all being binaries (Tsangarides et al., 2004; Lucatello et al., 
2005b). This binary mass transfer scenario is the same as 
that which is invoked to explain the Ba and CH stars (Iben 
& Renzini, 1983; McClure, 1984; McClure & Woodsworth, 
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1990; McClure, 1997). However, only about 1% of popu- 
lation I/II stars are Ba/CH stars, respectively (Tomkin 
et al., 1989; Luck & Bond, 1991). The CH stars are also 
carbon-rich but not as metal-poor as the CEMP stars, 
with [Fe/H] ~ — 1. The carbon-rich fraction of 1 per cent 
at higher metallicity is in stark contrast to the observed 
CEMP fraction of around 20 per cent. 

The mass-transfer scenario involves many processes that 
are not well understood. There are uncertainties associated 
with stellar evolution, particularly with respect to nucle- 
osynthesis in TPAGB stars. Carbon and s-process enhance- 
ments are thought to occur via third dredge up, but other 
processes may also play an important role in low-metallicity 
nucleosynthesis. These include hot-bottom burning (e.g. 
Iben, 1975; Boothroyd et al., 1993; Herwig, 2004), dual 
core flashes (also known as helium flash driven deep mix- 
ing) and dual shell flashes (helium flash driven deep mixing 
during a thermal pulse, Fujimoto, Iben, & Hollowell 1990; 
Schlattl, Salaris, Cassisi, & Weiss 2002; Cristallo, Straniero, 
Lederer, & Aringer 2007; Campbell & Lattanzio 2008), 
"extra mixing" on the first giant branch and perhaps on 
the TPAGB (e.g. Weiss, Denissenkov, & Charbonnel 2000; 
Nollett, Busso, & Wasserburg 2003; Eggleton, Dearborn, 
& Lattanzio 2008) and convective overshooting (Herwig, 
2000). 

The CEMP-formation scenario also requires knowledge 
of the physics of stellar interaction in binary systems. The 
primary star must transfer material to the secondary star 
which in turn might dilute and burn it. A wind mass- 
transfer scenario in wide binaries, e.g. by a mechanism 
similar to that of Bondi & Hoyle (1944), likely plays a 
role in CEMP formation. Closer binaries which undergo 
Roche-lobe overflow (RLOF) from a TPAGB star on to 
a less massive main-sequence star are expected to enter 
a common-envelope phase. Such stars would undergo few 
thermal pulses with little accretion on to the secondary 
(although see Ricker & Taam, 2008 for details of accretion 
in a common envelope and associated uncertainties). 

The fate of the accreted material is also uncertain. 
The molecular weight of accreted material is certainly 
greater than that of the secondary star, as carbon-enhanced 
TPAGB stars should also be helium rich. Accreted material 
should thus sink by the thermohaline instability (Stancliffe 
et al., 2007) but this may be inhibited by gravitational set- 
tling (Stancliffe & Glebbeek, 2008; Thompson et al., 2008). 
Furthermore, radiative levitation of some chemical species 
may be important (Richard et al., 2002a, b). When the sec- 
ondary ascends the first giant branch, its convection zone 
mixes any accreted material which may remain in the sur- 
face layers with material from deep inside the star. The 
surface abundance distribution depends on whether mate- 
rial has mixed deep into the star or not, because if it has 
it may have undergone some nuclear burning, the ashes of 
which are mixed to the surface. 

Two studies have considered population models in an at- 
tempt to reproduce the observed CEMP fraction of about 
20%. The models of Lucatello et al. (2005a) and Komiya 
et al. (2007) both concluded that in order to make enough 
CEMP stars the initial mass function (IMF) at low metal- 
licity must be significantly different to that observed in the 
solar neighbourhood. In particular, they enhanced the num- 
ber of intermediate-mass stars relative to low-mass stars - 
this has the effect of increasing the number of ~ 2 M stars 



which are responsible for the production of most of the car- 
bon. 

The Komiya et al. (2007) model differentiates between 
two metallicity regimes. In their model, stars with masses 
greater than 1.5 M Q undergo third dredge up irrespective of 
metallicity. For stars with [Fe/H] < —2.5 and mass less than 
1.5 M Q they invoke proton ingestion at the helium flash 
(the dual core flash) or at the first thermal pulse (the dual 
shell flash) as the source of carbon. This implies that the 
CEMP fraction should be smaller for [Fe/H] > -2.5, but 
the Stellar Abundances for Galactic Archeology (SAGA) 
database (Suda et al., 2008) shows that the CEMP fraction 
is approximately constant as a function of metallicity up to 
[Fe/H] « -2. 

A related problem is that of the nitrogen-enhanced 
metal-poor (NEMP) stars which have [N/Fe] > 0.5 and 
[C/N] < -0.5 as defined by Johnson et al. (2007). Such 
stars are expected to result from mass transfer in bina- 
ries with TPAGB primaries more massive than about 3 M Q 
in which hot bottom burning has converted most of the 
dredged-up carbon into nitrogen. The observed NEMP to 
EMP ratio is small, less than one in twenty-one according 
to Johnson et al. (2007) or less than 7% in the metallicity 
range [Fe/H] = —2.3 ± 0.5 according to the SAGA database 
(see Table 2). The Komiya et al. (2007) models, with an en- 
hanced number of intermediate-mass relative to low-mass 
stars, should make many more NEMP stars than are ob- 
served. 

The aim of this paper is to investigate which physi- 
cal scenarios are able to reproduce the CEMP and NEMP 
to EMP ratios, at metallicity [Fe/H] ~ —2.3, without al- 
tering the initial mass function. We combine a synthetic 
nucleosynthesis model with a binary population synthesis 
code to simulate populations of low-metallicity binaries (see 
Section 2). The power of the population synthesis approach 
is that it can efficiently explore the available parameter 
space. Much of the input physics is uncertain (as we have 
described above) but population synthesis allows us to ex- 
plore the consequences of these uncertainties by varying the 
model free parameters within reasonable bounds. We try to 
reproduce the observed CEMP and NEMP to EMP ratios, 
surface chemistry distributions, binary period distributions 
and chemical abundance correlations. It is thus a powerful 
tool to apply to this problem. 

In order to compare our models to observations we 
choose a subset of the SAGA database which corresponds 
to giants and turn-off stars as described in Section 3. The 
results of our simulations and comparison with the sam- 
ple of observations are given in Section 4. The implications 
of our results and outstanding problems are discussed in 
Section 5 while Section 6 concludes. 

2. Models 

In this section we describe our binary population synthesis 
model (Sections 2.1-2.3), initial distributions (Section 2.4), 
the choices of parameters for the various model sets 
(Section 2.5) and our criteria for selecting CEMP and 
NEMP stars (Section 2.6). 

2.1. Input physics 

Our binary population synthesis model is based on the syn- 
thetic nucleosynthesis models of Izzard et al. (2004) and 
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Izzard et al. (2006). Binary stellar evolution is followed ac- 
cording to the rapid binary stellar evolution (BSE) prescrip- 
tion of Hurley et al. (2002) in which detailed stellar evolu- 
tion model results are approximated by fitting functions. 
Coupled with a binary evolution algorithm which includes 
mass transfer due to both RLOF and winds, tidal circu- 
larisation and common envelope evolution, this approach 
allows the simulation of millions of binary stars in less than 
a day on a modern computer. 

Stellar evolution is augmented by a nucleosynthesis al- 
gorithm which follows the evolution of stars through the 
first, second and third dredge ups, altering surface abun- 
dances as necessary. This is mostly based on the Karakas 
et al. (2002) and Karakas & Lattanzio (2007, hereafter K02 
and K07 respectively) detailed models. 

We include a prescription for hot-bottom burning 
(HBB) in sufficiently massive AGB stars, M > 2.75 M 
at Z = 10~ 4 . The mass at which HBB switches on may 
be greater than our models suggest, e.g. > 5M in Weiss 
& Ferguson (2009, see their table B.4 for Z = 0.0005), 
3 - 4M Lau et al. (2009) or 3M in the latest mod- 
els by Karakas (2009, MNRAS submitted). These models 
use different input physics and/or nuclear reaction rates 
to the K02/K07 models on which our synthetic model is 
based. The impact of such changes on the number of CEMP 
stars is discussed in Section 4.1.1. Proton-capture reaction- 
rate uncertainties affect mainly hot-bottom burning stars, 
i.e. NEMP progenitors, rather than CEMP stars (see e.g. 
Izzard et al., 2007 for a discussion of the effect of proton- 
capture reaction rate uncertainties which affect Ne-Al in 
massive AGB stars). 

We model binary mass transfer by both stellar winds 
according to the Bondi-Hoyle prescription (Bondi & Hoyle, 
1944) and Roche-lobe overflow. Common-envelope evolu- 
tion follows the prescription of Hurley et al. (2002). We have 
updated some of the physical prescriptions in our model 
which are relevant to CEMP star formation. We describe 
below the most important changes to our binary code since 
Izzard et al. (2006). 

Our binary population synthesis model has been applied 
to a number of problems including the higher-metallicity 
equivalents of CEMP stars, the barium stars (Pols et al., 
2003) and CH stars (Izzard & Tout, 2004). Our model ap- 
proximately matches the observed Ba star to G/K giant 
ratio of ~ 1% (Luck & Bond, 1991). An extended version 
of our model was used by Bonacic Marinovic et al. (2008) 
to successfully model the eccentricities of the barium stars. 



2.1.1. Metallicity 

The Hurley et al. (2002) fitting formulae are limited to 
metallicities above and including Z = 10~ 4 and hence our 
stellar evolution model is not valid below this metallicity. 
Similarly, the K02/07 models extend down to Z = 10~ 4 
or, equivalently, [Fe/H] = —2.3. As a consequence we com- 
pare our models only to observations with [Fc/H] ~ —2.3 
(Section 3). We cannot compare our models to stars of sig- 
nificantly lower metallicity because our model lacks algo- 
rithms to describe phenomena such as proton ingestion at 
the first thermal pulse (see Section 2.1.6). 



2.1.2. First dredge up 

Abundance changes at first dredge up are interpolated from 
a grid of detailed stellar evolution models made with the 
STARS code. The STARS code was originally written by 
Eggleton (1971) and has been updated by many authors 
e.g. Pols et al. (1995) and Stancliffe & Eldridge (2009). 
The version used here employs the nucleosynthesis routines 
of Stancliffe et al. (2005), which follow forty isotopes from 
D to 32 S and important iron group elements. Model se- 
quences are evolved from the pre-main sequence to the tip 
of the red giant branch using 499 mesh points. Convective 
overshooting is employed via the prescription of Schroder 
et al. (1997) with an overshooting parameter of <5 v = 0.12. 
Thermohaline mixing on the RGB is included via the pre- 
scription of Kippenhahn et al. (1980). The diffusion coeffi- 
cient is multiplied by a factor of 100, following the work of 
Charbonnel & Zahn (2007). 

In single stars at low metallicity first dredge up has a 
small effect on surface abundances. However, in secondary 
stars which have been polluted by a companion first dredge 
up may either dilute accreted material which is sitting on 
the stellar surface or mix material from inside the star 
which has been been burned, depending on the efficiency of 
thermohaline mixing of the accreted material (Section 2.1.4; 
see also Stancliffe et al., 2007; Charbonnel & Zahn, 2007). 
We take such processes into account. A detailed description 
of our algorithm is given in Appendix A.l. 

2.1.3. Third dredge up 

Third dredge up is the primary mechanism by which car- 
bon made by helium burning is brought to the stellar sur- 
face in AGB stars. As stars evolve up the TPAGB their 
core mass M c increases and every Tip years a thermal pulse 
occurs. Once M c exceeds a threshold mass M Cim ; n third 
dredge up occurs with efficiency A, the ratio of the mass 
dredged up to the core growth during the previous inter- 
pulse phase. The values of A and M Cjm i n are fitted as a 
function of mass and metallicity to the detailed models of 
K02/K07. Without modification of this prescription single 
stars with initial mass greater than 1.25 M become carbon 
stars at a metallicity of Z = 10~ 4 . 

The correction factors AM Cim ; n and A m i n were intro- 
duced by Izzard et al. (2004) to enhance dredge up in low- 
mass stars relative to the detailed models 2 . They found 
that dredge up should occur earlier on the TPAGB and 
with greater efficiency than predicted by the K02 mod- 
els. With the parameter choices AM c ^ m ; n ~ — O.O7M 
and Amin ~ 0.8 — 37. 5Z the carbon-star luminosity func- 
tions in the Magellanic clouds are approximately fitted 
by the model. However, these parameters are poorly con- 
strained, especially at metallicities less than that of the 
Small Magellanic Cloud, so AM c ln i n and A m ; n should be 
considered free parameters at Z = 10~ 4 . 

We introduce a parameter, M env m i n , the minimum en- 
velope mass for third dredge up. This is 0.5 M by default, 
following solar- metallicity models (Straniero et al., 1997), 



2 These parameters modify the fits to the detailed models such 
that M Cjm i n — ► M c , m i n + AM Ciln i n and A — + max(A,A m i n ). A 
negative AM c , m in allows dredge up in lower initial-mass stars 
and a positive A m i n increases the amount of material dredged 
up once dredge up begins. 
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but we treat it as a free parameter. Recent detailed mod- 
els calculated with the stars code (Stancliffe & Glebbeek, 
2008) find third dredge up in an initially 0.9 M Q , Z = 10~ 4 
model which, when it reaches the AGB, has an envelope 
mass of only 0.31 M . Also, Stancliffe & Jeffery (2007) find 
third dredge up continues to occur even when the envelope 
mass drops below 0.5 M . In the latest models of Karakas 
(2009) dredge up is also found for small envelope mass 
(down to 0.1 M Q ) at Z = 10~ 4 . W e use this as justifica- 
tion of our decision to reduce M enVim ; n below 0.5 M Q but 
note that the value of Me nVj min is not accurately known. 

A choice of AM c _ m i n = —0.07, A m j n = 0.5 and 
M e nv,min = leads to dredge up in all stars which reach 
the TPAGB in the age of the galaxy, i.e. a minimum initial 
mass for third dredge up of ~ 0.8 M Q . 

In low-metallicity TPAGB stars dredge up of the 
hydrogen-burning shell is an important source of 1 C and 
14 N. We include an approximate prescription which well fits 
the K02/07 models (see Appendix A. 2). 



— Third dredge up parameters AM Cim i n = 0M Q , A m ; n = 
and Mcnv.min = 0.5 M Q which correspond to the 
detailed TPAGB models of K02/07. Enhanced third 
dredge up is simulated in some model sets by choosing 
a negative AM Cjln i n , positive A m ; n and zero M enVim i n . 

— Common envelope efficiency a = 1 according to the pre- 
scription of Hurley et al. (2002). The common-envelope 
structure parameter Ace is fitted to the detailed mod- 
els of Dewi & Tauris (2000). We do not include accre- 
tion on to the secondary star during the common en- 
velope phase by default but allow up to 0.05 M Q to be 
accreted in some model sets. We also consider the al- 
ternative common-envelope prescription of Nelemans & 
Tout (2005). 

— The 13 C pocket efficiency £13 is set to 1 by default as 
defined by Eq. A. 10 of Izzard et al. (2006). 

— The efficiency of the Companion Reinforced Attrition 
Process (CRAP, Tout & Eggleton, 1988) B is set to 
zero by default. 



2.1.4. Thermohaline mixing 

Most of our model sets assume thermohaline mixing of 
accreted material according to the prescription of Izzard 
et al. (2006) in which accreted material sinks and mixes in- 
stantaneously with the stellar envelope. The calculations of 
Stancliffe et al. (2007) suggest this is reasonable in some 
cases. However, gravitational settling prior to accretion 
may prevent thermohaline mixing (Thompson et al., 2008; 
Bisterzo et al., 2008). In order to account for both possibil- 
ities, either efficient and instantaneous or highly inefficient 
thermohaline mixing, we have also run models in which the 
accreted material remains on the stellar surface (regard- 
less of its molecular weight) until mixed in by convection. 
Recent calculations show that the situation is somewhat 
more complicated than either of the extremes we test here 
(Stancliffe & Glebbeek, 2008). 

2.1.5. Parameter choices 

Our binary nucleosynthesis model has many free parame- 
ters, some of which have been constrained by previous stud- 
ies, some which have not. We list here the most important 
parameters and our default choices. 

— Abundances are solar-scaled with a mixture according 
to Anders & Grevesse (1989). We do not include an a- 
element enhancement. Most of our models have a metal- 
licity of Z = 10~ 4 (equivalent to [Fe/H] = -2.3). 

— Wind mass-loss rates are parameterised according to the 
Reimers formula with r] — 0.5 on the first giant branch 
and Vassiliadis & Wood (1993) on the AGB (as in K02). 
We modulate the AGB mass-loss rate with a factor /vw 
which is one by default. We apply a correction APyw, 
zero by default, to the Mira period relation used in the 
Vassiliadis & Wood prescription to simulate a delayed 
superwind on the AGB. We also consider both Reimers 
and Van Loon mass- loss rates on the AGB (Reimers, 
1975; van Loon et al., 2005). Appendix B describes the 
mass-loss formulae in detail. 

— The Bondi-Hoyle accretion efficiency factor ckbh = 3/2. 
We also consider «bh = 5, however unphysical this may 
be, to simulate enhanced stellar-wind mass transfer. 



2.1.6. Missing physics 

Our synthetic models do not include any extra mixing 
which may be responsible for the conversion of 12 C to 
13 C and 14 N in low-mass stars (see e.g. Nollett, Busso, & 
Wasserburg, 2003; Busso et al., 2007; Eggleton, Dearborn, 
& Lattanzio 2008 and references therein). Also, we do not 
include any prescription which describes mixing events in- 
duced by proton ingestion at the helium flash (the dual core 
flash) or during thermal pulses (the dual shell flashes, both 
of which are also known as helium-flash-driven deep mix- 
ing; Fujimoto et al., 1990; Weiss et al., 2004; Cristallo et al., 
2007; Campbell & Lattanzio, 2008). Current stellar models 
suggest these events occur only if [Fe/H] < —3 while our 
models have [Fe/H] = —2.3 and our observational sample 
includes stars with [Fe/H] = — 2.3±0.5. The latest results of 
the Teramo group show that at the lowest masses, around 
0.8 M Q for [Fe/H] = -2.3 (with Z Q ~ 0.01), a dual shell 
flash may occur at the beginning of the TPAGB with sig- 
nificant C, N and s-process production (Cristallo, private 
communication) - we leave the analysis of such a case to 
future work. 



2.2. An example CEMP system 

We preempt the results of Section 4 with an example of our 
synthetic stellar evolution algorithm for a CEMP system as 
shown in Fig. 1. Initially it has Mi — 2 M Q , M2 = 0.75 M Q , 
P = 3700 days (a = 1400 R Q ), Z = 10~ 4 , e = and oth- 
erwise has the default physics as described above. The pri- 
mary mass is chosen to illustrate the case of maximum car- 
bon accretion. After 880 Myr the primary evolves onto the 
TPAGB. Third dredge up increases its surface carbon abun- 
dance to [C/Fe] = +3.2 (Fig. 16). Of its wind 0.11 M Q is 
accreted by and mixed into the main-sequence secondary 
such that its mass increases to 0.86 M Q (Fig. la), it has a 
carbon abundance [C/Fe] = +2.3 and nitrogen abundance 
[N/Fe] = +0.35 (Fig. 16). The binary orbit expands to 
a period of 5500 days (Fig. la). At this moment the sec- 
ondary becomes a CEMP star, although in our population 
synthesis we do not count it as such until it has evolved to 
logg < 4 (see Section 2.6), which occurs after 10.59 Gyr. 
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Fig. 1. An example CEMP system with our default parameter choices as in model set A. Initially Mi = 2M0, M2 = 0.75 M© and 
P = 3700 days. Panels a and b show the evolution of the primary up the AGB and accretion onto the secondary. Panels c and d 
show the evolution of the secondary through the CEMP phase. In panels a and c the left ordinate denotes the mass of the stars in 
Mq (solid and dashed lines for the primary and secondary respectively), the right ordinate gives the orbital period in days (dotted 
lines). In panel b the left ordinate shows [C/Fe] while in panel d the left axis shows [N/Fc]. In all panels the abscissa shows the 
age of the system. Labels in italic show the evolution phases, labels in Roman type show the dredge up, mixing and accretion 
phases, and arrows in panels c and d indicate when the star is tagged as a CEMP star according to our selection criteria (older 
than lOGyr, log g < 4 and [C/Fe] > 1, see Section 2.6). 



Subsequently the secondary also ascends the giant 
branch twice (Fig. lc and d). During the first ascent, dredge 
up reduces the surface abundance of carbon only slightly 
but increases [N/Fc] to +1.2 when burnt material is mixed 
to the surface. The orbit expands when the secondary loses 
0.15 Mq through its stellar wind at the tip of the giant 
branch (Fig. lc). At 11.3 Gyr it too ascends the AGB with 
a mass of 0.7 M Q and ceases to be a CEMP star when its 
envelope is lost in a wind and it becomes a white dwarf. 
The system is then a pair of carbon-oxygen white dwarfs 
with a period of about 20 years. 

Further examples of our stellar evolution and nucleosyn- 
thesis algorithms and comparisons to detailed stellar evo- 
lution models are given in Izzard et al. (2004, 2006). 



2.3. Population synthesis 

Each of our population synthesis simulations consists of N 3 
stars in lnM 1 -lnM2-lna parameter space, where M lj2 are 
the initial masses of the primary and secondary, a is the 
initial separation and N = 128. We then count the number 



of stars of a particular type according to the sum 

Mi, max M 2 , max a max t max 

nt yP c=Sj2 E J2Y,S(typc)^SM 1 8M 2 SaSt, (1) 

Mi >min M 2imin a min t min 

where S is the star formation rate, ^ is the initial dis- 
tribution function and (5(type) = 1 when a star is of the 
required type and zero otherwise. The grid cell size is given 
by 5Mi -5M2 -Sa while the timestep is 5t. We further assume 
that \1/ is separable, 

* = V(MiMM 2 ) X (a), (2) 

where the functions ip(Mi), and x( a ) are the initial 

distributions of Mi, M2 and a respectively (see Section 2.4). 
The star formation rate S is a function of time or, given an 
age-metallicity relation, metallicity. Because we calculate 
ratios of numbers of stars with the same metallicity, i.e. 
the same age and same S, the star formation rate simply 
cancels out. 

We set the limits of our population synthesis grid as 
follows: 

— Mi :m i n = 0.7 M Q - stars less massive than this do not 
evolve off the main sequence within the lifetime of the 
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Universe. Mi jmax = 8M - we do not consider massive 
stars which end their lives as supernovae. 

— -^2,min = 0.1 Mq - There is no simple limit on M2 >m i n 
as a secondary star may accrete an arbitrary fraction 
of the mass lost from its companion. However, with 
the Bondi-Hoyle accretion formalism the accreted mass 
is typically less than 0.2 M©, so we catch all possible 
CEMP stars (which must have masses around 0.8 M©). 
-^2,max = 0.9 M© - stars more massive than this have 
already evolved to white dwarfs after ~ lOGyr and so 
cannot be CEMP stars. 

— a m i n = 3R©, typically a max = 10 5 R©. The upper limit 
is chosen to include all CEMP stars made in our models. 
Our assumption is that all stars are made in this range 
of separations, i.e. a binary fraction of 100%. In reality 
some systems are wider and/or single. This can easily 
be accommodated by lowering the binary fraction. 

— t mm = lOGyr - Our stars must be old halo stars, so 
a lOGyr limit is reasonable. Our results are not sensi- 
tive to changes of ±2Gyr in this limit. The age of the 
Universe gives the upper limit i max = 13.7 Gyr. 

Selection criteria for £(type) are given in Section 2.6. 

2.4. Stellar distributions 

The choice of distributions of initial primary mass Mi , ini- 
tial secondary mass Mi (or, alternatively, q = M2/M1) 
and initial separation a (or initial period P) affects the fi- 
nal number counts and distribution of CEMP parameters. 
Unfortunately, in the Galactic halo all the relevant initial 
distributions are unknown. We are forced to assume solar 
neighbourhood distributions: 

— The primary mass distribution ip(M\) is the initial mass 
function of Kroupa, Tout, & Gilmore (1993, KTG93). 

— The secondary mass distribution (^{M?) is flat in q = 
M 2 /Mi, i.e. any mass ratio is equally likely. 

— The separation distribution x{ a ) ls flat in In a between 
a m in and a max (i.e. x( a ) ~ !/«)• 

— We start all binaries in circular orbits, i.e. eccentricity 
e = 0, and assume a binary fraction of 100%. 

2.5. Model sets 

As described in Section 2.1.5, many physical parameters as- 
sociated with binary evolution are uncertain. Table 1 lists 
the most important parameter sets we consider and how 
they differ from our default model set (set A, with param- 
eter values defined in Section 2.1.5, see also Table D.l for a 
list of all the model sets considered, including those not dis- 
cussed in detail). In some model sets, e.g. Ap5 and Ap7, we 
altered the CEMP selection criteria rather than the physical 
parameters. All our model sets use the initial distributions 
outlined in Section 2.4 (except model set Ae5 with initial 
eccentricity e = 0.5). 

2.6. Model selection criteria 

Stars are selected from our model population as metal poor 
(EMP) giants if their surface gravity log 10 g < 4 and they 
are older than 10 Gyr corresponding to the approximate age 
of the Galactic halo. We then define the following subtypes: 



IVlULlcl ScL 


I 11 V &ll_.dl pdl dlllc Lcl 

(differences from model set A) 


A 

Ap5, Ap7 


[0/J H eJ min = U.5 and 0.7 respectively 


A 1 A Hi 

Al, Az 


£13 =0.1 and 0.01 respectively 


B 


Bondi-Hoyle efficiency «bh = 5 


c 


common envelope accretion 0.05 Mq 


D 


no thermohaline mixing 


E 


AM c , min = -0.07 M Q , A min = 0.8 


F 


-Menv.min = Mq 


G 


AM c , min = -0.1 Mq, A min = 0.8, 

M cnv , min = 0M Q , £i 3 = 0.1 


H 


as model G, with: 
common envelope accretion 0.05 Mq 
no thermohaline mixing 



Table 1. Physical parameters corresponding to the most im- 
portant of our binary population models. A full list of all our 
model sets can be found in Table D.l. 



CEMP Surface carbon abundance [C/Fe] > [C/Fe] mi „, 
where [C/Fe] m i n = 1.0 by default. If in addition 
[Ba/Fc] > 0.5 the star is classified as CEMP-s. 

NEMP [N/Fe] > 1.0 and [C/N] < -0.5. The former is 
somewhat more restrictive than the criterion [N/Fe] > 
0.5 applied by Johnson et al. (2007) and more consis- 
tent with the CEMP criterion. We find that with the 
Johnson et al. (2007) criterion normal giants that have 
undergone strong CN cycling, presumably due to extra 
mixing on the giant branch which is not included in our 
models, enter our observational selection (see Section 3). 
Therefore we define the NEMP class to contain only 
stars affected by third dredge-up and HBB and subse- 
quent mass transfer, i.e. the nitrogen-rich equivalents of 
CEMP stars. 

FEMP Surface fluorine abundance [F/Fe] > 1.0 (Lugaro 
et al., 2008). 

These three subtypes are not mutually exclusive. It turns 
out that in our models the FEMP and CEMP sub- 
types nearly coincide (see Section 4.1.3). The CEMP and 
NEMP classes also partially overlap, these are designated 
as CNEMP. 

3. Observational database and selection criteria 

Our database of observed EMP stars is based on 2376 ob- 
servations of 1300 stars in the SAGA database 3 (Suda et al., 
2008). When observed values of a parameter are available 
from multiple sources for one star, we simply use the arith- 
metic mean of these values in our database (see Appendix E 
for details). We ignore observations which provide only an 
upper or lower limit. We select stars which correspond to 
our model EMP giants and turn-off stars as follows: 

1. The observed star must have metallicity in the range 
[Fe/H] = -2.3 ± 0.5 dex. Fig. 2a shows that the num- 
ber distribution of stars in this range varies by a factor 
of two with no clear trend in number as a function of 
[Fe/H] . 

2. The star must be a giant or sub-giant, i.e. log g < 4.0. 
Fig. 2b shows the distribution of the number of stars as 
a function of log g. 

3 With the addition of log g values from Lucatello et al. (2006, 
Beers, private communication). 
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Fig. 2. Selected distributions from our sample of the SAGA database. All stars in our selection have log g < 4 and [Fe/H] = 
—2.3 ± 0.5. panels a and b show the distributions of [Fe/H] (the metallicity distribution function) and \ogg respectively for the 
sets of 1) all stars in our selection and 2) those stars in our selection with a carbon measurement. The horizontal error bars, and 
hence the bin widths, are typical la errors for a single observation, ±0.15 dex in [Fe/H] and ±0.25 dex in logg. 



This selection leaves us with 373 stars, of which 308 have 
a measured carbon abundance and 96 have [C/Fe] > 1 
these are our CEMP stars 4 . The CEMP to EMP fraction, 
for stars with measured carbon, is thus 31%. Fig. 2a shows 
that whether a star has measured carbon, or not, does not 
depend on metallicity but is sensitive to logg. In particular, 
90% of stars with logg < 3.5 have carbon measurements, 
whereas about half of the highest-gravity stars (logg > 
3.5) do not. Presumably this is because of observational 
difficulties, these stars being relatively dim and hot. 

If we assume that stars which have no carbon measure- 
ment actually have no carbon enhancement, i.e. [C/Fe] < 1, 
the CEMP/EMP ratio drops to 26%. The latter statistic is 
best for comparison with our models, but it rests on a possi- 
bly dubious assumption: are stars with no carbon measure- 
ment really not enhanced in carbon? Certainly the CEMP 
fraction depends on how it is counted, as shown in Table 2. 
The CEMP fraction also varies depending on the survey un- 
der consideration: Frebel et al. (2006), Cohen et al. (2005) 
and Lucatello et al. (2006) find 9, 14 and 21% respectively. 

Application of the NEMP criteria 5 , [N/Fe] > 1 and 
[C/N] < -0.5, leaves us with zero NEMP stars in the 
[Fe/H] range considered. This is to some extent a statis- 
tical fluke, because NEMP stars are found both at higher 
(in the case of HE0400-2040) and especially at lower metal- 
licity. However, even when counting all stars regardless of 
[Fe/H] the number of NEMP stars is small compared to 
the total number of EMP or CEMP stars (see Table 2). 
Use of the less restrictive [N/Fe] > 0.5 criterion would in- 
troduce seven NEMP stars into the sample, but these are 
mostly carbon-depleted giants that have presumably un- 
dergone strong CN-cycling during their RGB evolution, as 
discussed in Section 2.6. 

Our selected sample may be biased, as shown in Fig. 
3b. The CEMP/EMP ratio increases slightly as logg de- 
creases but a constant value is consistent with the error 
bars. The CEMP fraction may increase slightly with metal- 



4 When [C/Fe] is not available but [CH/Fe] is, we use [CH/Fe] 
as a proxy for [C/Fe]. 

6 We use [CN/Fe] or [NH/Fe] as proxies for [N/Fe] when 
[N/Fe] is not available in the SAGA database. 



licity (Fig. 3a) at least over the narrow range we consider 
([Fe/H] = -2.3±0.5). If dual-core and/or dual-shell flashes 
occur exclusively for [Fc/H] < — 2.5 and are responsible for 
the formation of most CEMP stars we expect the CEMP 
fraction to drop as [Fe/H] increases beyond —2.5. This is 
not seen. We note that the SAGA database was not de- 
signed to be complete in any statistical sense. We await a 
more complete census of metal-poor stars before definite 
conclusions on the CEMP fraction and its dependence on 
metallicity and evolutionary stage, can be drawn. 

4. Results and comparison with observations 

The ratios of the number of CEMP and NEMP stars to 
EMP stars as selected from our model sets of Table 1 are 
shown in Table 3. Our results fall into two categories: 

1. Most of our model sets have a CEMP to EMP ratio of 
2 — 4%. This is not much larger than the 1% of more 
metal-rich giants which are CH stars. We describe below 
one of these, our default model set (model set ^4), in 
detail. 

2. A few sets come close to reproducing the observed 
CEMP to EMP ratio. These are the sets with some com- 
bination of the following: (enhanced) third dredge up in 
low-mass stars, no thermohaline mixing and accretion 
in the common envelope phase. They represent a corner 
of the parameter space which may be considered rather 
extreme, though not unfeasible. We describe sets G and 
H in detail below. 

4.1. Model set A: Default physics 

Our default model set - set A - represents a choice of phys- 
ical parameters which could be described as conservative. 
The parameters are not controversial or extreme. As such 
they are a good starting point for our analysis of the CEMP 
problem. 

First of all, the CEMP to EMP ratio in our default 
model set is 2.3%. This is clearly at odds with the observed 
9 — 25% ratio, especially if one factors in a binary fraction 
smaller than unity. However, the number of NEMP stars 
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Metallicity 


logs 


EMP 


CEMP 


CNEMP 


NEMP 




CEMP/EMP 


Additional Selection Criteria 


[Fe/H] 










[N/Fe] > 1 


[N/Fe] > 0.5 


fraction 




-2.3 ±0.5 


< 4 


308 


96 








31% 


only stars with C measured 


-2.3 ±0.5 


< 4 


373 


96 








26% 


all EMP stars 


-2.3 ±0.5 


< 4 


104 


88 








7 




only stars with C and N measured 


-2.3 ±0.5 


< 4 


373 


96 








7 


26% 


f C and N measured for (C)NEMPs, 
y C measured for CEMPs, 


all 


< 4 


77Q 




A 


1 A 


99 
zz 


1(5/0 




-2.3 ±0.5 


all 


479 


115 








7 


24% 




all 


all 


1366 


177 


6 


17 


23 


13% 




< 2.0 

< 2.0 


giants 


132 

270 


11 

58 








9 ±2% 
> 21 ± 2% 


Frebel et al. (2006) 
Lucatello et al. (2006) 


< 2.0 


giants 












14 ± 4% 


Cohen et al. (2005) 



Table 2. Observed numbers of EMP, CEMP, CNEMP and NEMP stars from the SAGA database for our various selection criteria. 
The first line of data counts only stars which have a [C/Fe] measurement. The second line of data additionally counts stars without 
measured [C/Fe] as EMPs, i.e. it is assumed that [C/Fe] < 1 for these stars. In the next four lines stars with both measured C 
and N are used to count (C)NEMP stars while stars with C measured count the CEMP stars. 

Also shown are estimates of the CEMP/EMP ratio from Cohen et al. (2005), Frebel et al. (2006) and Lucatello et al. (2006). 



0.5 |— r 




-2.8 -2.6 -2.4 -2.2 -2 -1.8 1 2 3 4 



[Fe/H] log g 

Fig. 3. The CEMP/EMP ratio for our selection from the SAGA database ([Fe/H] = -2.3 ± 0.5 and log 3 < 4) as 

a) a function of metallicity [Fe/H] and 

b) a function of gravity log g. 

Error bars are based on Poisson (y/N) statistics and bin widths are as in Fig. 2. 



Model Set 


CEMP/EMP % 


NEMP/EMP % 


FEMP/EMP % 


CEMP-s/CEMP % 


A 


2.30 ± 0.04 


0.267 ± 0.007 


2.03 


28.3 


Ap7 


2.99 ±0.04 


0.267 ± 0.005 


2.03 


21.8 


Ap5 


3.46 ± 0.04 


0.267 ± 0.005 


2.03 


18.8 


Al 


2.30 ± 0.04 


0.267 ± 0.007 


2.03 


86.3 


A2 


2.30 ± 0.04 


0.267 ± 0.007 


2.03 


94.1 


B 


3.05 ± 0.04 


0.323 ± 0.006 


2.74 


47.2 


C 


2.94 ±0.04 


0.311 ±0.007 


2.36 


29.1 


D 


4.21 ±0.04 


0.409 ± 0.004 


3.74 


63.1 


E 


2.90 ± 0.04 


0.267 ± 0.007 


2.43 


43.5 


F 


6.47 ±0.03 


0.267 ± 0.006 


5.16 


22.3 


G 


9.43 ± 0.04 


0.266 ± 0.006 


7.81 


94.1 


H 


15.52 ±0.07 


0.426 ± 0.004 


13.5 


97.2 



Table 3. Percentage of CEMP, NEMP and FEMP (sub-)giants relative to total EMP giants in our model binary populations of 
Table 1 (see Section 2.6 for selection criteria and Table D.2 for the full set of results). The errors convey Poisson statistics only. 
The final column gives the number ratio of CEMP-s to CEMP stars. 



in this simulation is small (0.3% of EMP stars), which does 
agree with the observations. In this section we examine var- 
ious properties of our default population with a view to 
later sections which improve the match between our mod- 
elled CEMP to EMP ratio and the observations. 



4.1.1. CEMP initial parameter space 

The regions of the initial Mi — M<i — a or P parameter space 
which form CEMP stars are shown in Figure 4. The distri- 
butions are time- weighted, as in Eq. 1. The vast majority 
of our CEMP stars form via the wind-accretion channel 
with a typical Mi ~ 1.2 - 1.5 M Q , M 2 ~ 0.8 M , sep- 
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Initial M l /M i 



0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 
Initial M,/M, 



2 /lVlg 




3.2 3.4 3.6 3.8 
log Initial Separation/R,., 



4.2 




3.5 4 4.5 

log Initial Period/days 



Fig. 4. Time- weighted distributions of initial masses, Mi, M2, initial separation a and initial period P in binaries which lead to the 
formation of a CEMP star in our default model set A. In this model set almost all our CEMP stars are formed by accretion from 
primaries with M > 1.25 Mq. The vertical scale, which counts the number of stars per bin, is linear and arbitrarily normalized to 
peak at one. 



aration around 10 3 — 10 4 R0 or, equivalently, an orbital 
period of 10 3 — 10 5 days. The initial mass Mi is limited 
at the low end by the minimum mass for third dredge up 
(~ 1.2 M©) and the distribution peters out at the high end 
mainly as a result of the IMF. Stars with Mi > 2.7 M Q 
undergo hot-bottom burning which results in [C/N] < — 1. 
These also classify as (C)NEMP stars. As mentioned in 
Section 2.1 some uncertainty exists regarding the mass of 
HBB onset which may be greater than 2.7 M©. Because of 
the rapid drop in the IMF with increasing mass, as shown 
in Figure 4a, the number of CEMP stars affected by this 
uncertainty is small. An increase in the HBB-onset mass 
reduces the number of NEMP stars. 

The secondary mass, M 2 ~ 0.8 M©, is the mass expected 
for a star that is approximately 10 Gyr old. The shortest pe- 
riod (or separation) binary which forms a CEMP is limited 
by the Roche limit. Closer binaries pass through a common- 
envelope stage with little accretion on the secondary. The 
efficiency of wind mass transfer drops as the initial separa- 
tion increases. Beyond about 10 4 R© the secondary accretes 
too little carbon to become a CEMP. 

4.1.2. The distribution of \ogg and potential selection effects 

In Fig. 5 we compare the distribution of log g in model set A 
CEMP stars to our selection from the SAGA database (see 
Fig. 3b). It is hard to understand why the distributions dif- 





Lj|| 


(-1 ■ 









12 3 4 



Fig. 5. The distribution of logg in CEMP stars taken from 
our default stellar population model A (filled histogram) vs our 
CEMP selection from the SAGA database (open histogram with 
Poisson error bars). 

Note that in this and the similar plots that follow the area un- 
der the graph (which is the total number of stars) is normalized 
such that it is the same for both the observations and our model 
stars. 

fer without invoking selection effects such that low-gravity 
CEMP giants are preferred (see Section 3). We have made 
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1 1.5 2 
[ 14 N/ 56 Fe] 



Fig. 6. The distribution of [C/Fe] in our default CEMP popu- 
lation (model set A, filled histogram) compared to observations 
(open histogram with Poisson error bars). 



Fig. 7. The distribution of [N/Fe] in our default CEMP popu- 
lation A (filled histogram) compared to observations (open his- 
togram with Poisson error bars). 



no attempt to model such selection effects but we do not 
expect this to strongly affect our model predictions. The 
log g distribution of non-CEMP stars in model set A has 
a nearly identical shape to the CEMP stars, so that the 
CEMP/EMP ratio for our default population is indepen- 
dent of log g and hence the overall CEMP fraction pre- 
sented in Table 3 is quite robust. We also find that, at least 
in our default model as well as other models that include 
thermohaline mixing of accreted material, there is little sur- 
face abundance evolution in our CEMP stars after accretion 
(with the exception of nitrogen, which increases somewhat 
at first dredge up) so that the abundance distributions we 
present in the following subsections are also hardly affected 
by this selection bias. 

Finally, we note that the peak in the model distribution 
at logg ~ 3.2 is due to horizontal branch stars, most of 
which have effective temperatures close to 10,000 K in our 
models. The SAGA database does not contain horizontal 
branch stars hotter than ~ 7,000 K probably because hotter 
stars are selected against. If we select only red horizontal 
branch stars cooler than 7,000 K a small peak at logg ~ 2 
remains. This may correspond to the tentative peak seen 
in the observed distribution. 
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Fig. 8. The distribution of [N/Fe] versus [C/Fe] in our default 
CEMP population, model set A (darker grey indicates a larger 
density of stars). The vertical dashed line indicates our CEMP 
selection criterion ([C/Fe] > 1) and the diagonal dashed line 
shows our NEMP selection criteria ([N/Fe] > 1 and [N/C] > 
0.5). Observed CEMP stars are indicated by crosses. 



4.1.3. Carbon, nitrogen and fluorine 

The distribution of carbon in the CEMP population of 
model set A matches reasonably well the observed distribu- 
tion, as shown in Fig. 6. The observed distribution shows 
somewhat higher carbon enhancements on average, while 
our models predict more stars in the range 1 < [C/Fe] < 1.5 
and none above [C/Fe] = 2.5, compared to about ten ob- 
served (11 per cent of the sample). This may be partly ex- 
plained by our assumption of complete thermohaline mix- 
ing: we assume the entire star mixes, when in reality only 
a fraction mixes. 

The picture becomes much less favourable when we com- 
pare the distribution of nitrogen, as shown in Figs. 7 and 8. 
Our default population contains few (C)NEMPs (above the 
dashed line in Fig. 8; note that NEMPs with [C/Fe] < 1.0 
are not shown in this figure), which agrees with both our 
observational sample and that of Johnson et al. (2007). 
However, most observed CEMP stars are enhanced in nitro- 



gen by 1 — 1.5 dex, which in fact coincides with a dearth of 
model CEMP stars. Our model includes three mechanisms 
for increasing the nitrogen abundance: first dredge up, hot- 
bottom burning and third dredge up of the hydrogen burn- 
ing shell. The latter is only a small effect and while first 
dredge up enhances [N/Fe] by typically 0.5 dex, this cannot 
reproduce the 1 dex nitrogen enhancements seen in our ob- 
servational sample. On the other hand, hot-bottom burning 
converts most of the dredged-up carbon into nitrogen and 
thus results in much larger nitrogen enhancements. If HBB 
were more effective than we assume it would only raise the 
number of (C)NEMP stars, in contradiction with the obser- 
vations. We must therefore assume that either some kind 
of extra- mixing mechanism or a dual core/shell flash is re- 
sponsible. This is beyond the scope of our present model. 

Another indication of a missing ingredient in our models 
comes from the carbon isotopic ratio. The 12 C/ 13 C ratio in 
our models is always large (10 2 — 10 4 ), whereas the observed 
ratio is generally less than the solar ratio (around 90) , from 
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Fig. 9. The distribution of [F/Fe] versus [C/Fe] in our model 
set A. The horizontal dashed line indicates our FEMP selection 
criterion ([F/Fe] > 1) and the vertical dashed line indicates 
our CEMP selection criterion ([C/Fe] > 1). The position of the 
fluorine-rich CEMP star HE 1305+0132 is also plotted. 
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Fig. 10. Sodium distribution in our default model set A (filled 

histogram) compared to observations (open histogram with 
Poisson error bars). 



the equilibrium value of four up to about 50 (Ryan et al., 
2005). Again, this may result from extra mixing or a dual 
shell/core flash. 

Fluorine was recently measured in the CEMP star HE 
1305+0132 (Schuler et al., 2007) and the halo planetary 
nebula BoBn 1 (Otsuka et al., 2008). Lugaro et al. (2008) 
pointed out that fluorine is made in the progenitors of 
CEMP stars and therefore that most CEMP stars should 
be FEMP stars (EMP stars with [F/Fe] > 1). Fig. 9 con- 
firms this for model set A. Our model struggles to reproduce 
the [F/Fe] ratio of HE 1305+0132 although this star may 
represent the tail of the fluorine distribution. Future obser- 
vations of fluorine in CEMP stars should reveal this pop- 
ulation. Also, the abundance may be overestimated (Abia 
et al., 2009). 

4.1.4. Sodium and other light elements 

Sodium is enhanced at the surface of TPAGB stars through 
both third dredge up and hot-bottom burning. In our model 



set A most stars dredge up little sodium, leading to the 
peak at [Na/Fc] ~ +0 in Fig. 10. A few of our stars, 
with relatively massive AGB primaries, dredge up enough 
sodium that the secondary reaches up to [Na/Fc] ~ 1.5. 
The most massive AGB primaries undergo HBB and give 
rise to the small peak seen at [Na/Fe] ~ 2 — 2.5: these stars 
are CNEMP stars. Our model at least matches the range of 
observed stars, although sodium-enriched stars are greatly 
underrepresented. Note that there is great uncertainty in 
the yield of sodium from massive AGB stars, and hence in 
our model CNEMP sodium abundances, because of reac- 
tion rate uncertainties (Izzard et al., 2007). 

Observations of CEMP stars show an overabundance of 
oxygen and other a-elements typical of the halo, i.e. about 
+0.4 dex. Our models include no explicit a-enhancement 
but do show enhancements in oxygen and magnesium be- 
cause of third dredge up (up to +0.5 dex and +0.6 dex re- 
spectively) . The few observations of oxygen in CEMP stars 
range from to +2 dex while magnesium is enhanced by up 
to +1.5 dex. It is clear that our models struggle to repro- 
duce these stars, especially because they are often giants 
which should be well mixed. 

There have been several recent lithium abundance mea- 
surements in CEMP stars (see e.g. Thompson et al., 2008 
and Roederer et al., 2008). Modelling lithium nucleosyn- 
thesis is too complicated for our synthetic models. Detailed 
discussions can be found in Roederer et al. (2008) and 
Stancliffe (2009). 

4.1.5. The heavy elements 

Model A yields a CEMP-s/CEMP ratio of only 28 per 
cent, smaller than is observed: of the 47 CEMP stars in 
our SAGA database selection with both carbon and bar- 
ium measured 6 44 have [Ba/Fe] > 0.5 (94%), while Aoki 
et al. (2007) report an s-rich fraction of 80%. All our model 
CEMP stars are enriched in s-process elements, but most 
stars in this model have < [Ba/Fe] < 0.5 and are thus 
not classified as CEMP-s. This is because the assumed 13 C 
pocket efficiency £13 = 1 gives such a high neutron expo- 
sure that the s-process distribution is pushed to the lead 
peak (Gallino et al., 1998). A decreased £13 (model sets Al 
and A 2) gives larger barium abundances and hence a larger 
CEMP-s fraction (see Table 3). 

The need to decrease the 13 C efficiency at low metallic- 
ity, to £13 ~ 0.1 for [Fe/H] < —1, was shown by Bonacic 
Marinovic et al. (2007) on the basis of [Pb/hs] ratios in lead 
stars. A comparison between model set A and the observed 
heavy element abundance distribution is shown in Fig. 11. 
We find a best match to the [Pb/Ba] ratio for model Al 
with £13 = 0.1. Although it is unsatisfactory that both the 
CEMP-s to CEMP ratio and the s-abundance ratio distri- 
butions depend quite sensitively on a free parameter in the 
model, we find a reasonable match to all these constraints 
for a single value of £13. 



4.1.6. Orbital periods 

Comparison of observed CEMP orbital periods with our 
models is difficult for two reasons. First, the number of stars 



We could assume that stars without measured barium have 
little barium, in which case the number of CEMP stars is 96 and 
the s-rich fraction is 44/96 = 46%. 
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Fig. 11. Heavy elements in CEMP stars from model sets A, Al and A2 with the 13 C efficiency parameter £13 = 1, 0.1 and 0.01 
respectively. Our observational sample is shown by the filled histograms (with Poisson error bars) while the model sets A, Al and 
A2 are plotted with long-dashed, short-dashed and dot-dashed lines respectively. The distributions are all normalized to peak at 
one. We choose Y, Ba and Pb to represent the s-process peaks. 



Object 


Period 


e 


Reference 


CS22948-027 


505d 


0.3 


Preston & Sneden (2001) 


CS22948-027 


426.5d 


0.02 


Barbuy et al. (2005) 


HD5223 


755.2d 





McClure & Woodsworth (1990) 


HD224959 


1273d 


0.179 


McClure & Woodsworth (1990) 


CS22956-028 


1290d 


0.22 


Sneden et al. (2003) 


CS22942-019 


2800d 


0.19 


Preston & Sneden (2001) 


LP 625-44 


12y 




Aoki et al. (2000) 



Table 4. Periods and eccentricities of CEMP stars in binaries 
from our selection ([Fc/H] = —2.3 ± 0.5, log 10 g < 4). Note that 
two solutions are available for CS22948-027. 



with known periods is small despite the fact that many 
are binaries - from our selection there are only six stars 
with orbital solutions (see Table 4). Second, long periods 
are difficult to measure. Any period longer than about ten 
years, which is approximately the time for which surveys 
have been ongoing, is likely to remain unmeasured for some 
time to come. With this in mind, Figure 12 shows the dis- 
tribution of periods from model set A compared to the six 
observed stars. 

The lower limit to the model period distribution is set 
by systems that are just wide enough for the AGB pri- 
mary to avoid filling its Roche lobe when the star has 
its maximum radius. The shortest-period CEMP binaries, 
CS22948-027 and HD5223, are not compatible with our 
model set A. Their short periods and small eccentricities (if 
we adopt the Barbuy et al., 2005 result regarding CS22948- 
027) suggest they underwent tidal circularisation, and per- 
haps RLOF and common-envelope evolution. This does not 
explain their carbon enhancement. The problem is remi- 
niscent of the barium stars, which have too large eccen- 
tricities and too short periods compared to models (Pols 
et al., 2003). Various explanations have been proposed to 
deal with this problem which build on the fact that AGB 
stars have very extended, outflowing atmospheres such that 
the canonical distinction between RLOF (when the stellar 
surface fills the Roche lobe) and wind accretion (when it 
does not) becomes blurred (e.g. see Frankowski & Jorissen, 
2007, Bonacic Marinovic et al., 2008). Our understanding 
of binary evolution in this transition region is still poor. An 




2 3 4 5 6 7 
log 10 (Period/days) 



Fig. 12. Period distribution in our CEMP model sets A, G and 
H (solid lines, dashes and dot-dashes respectively) vs CEMP 
period measurements from our selection of the SAGA database 
(crosses) . 



alternative solution may be accretion during the common- 
envelope phase (Ricker & Taam, 2008), which we explore 
in model H (see Section 4.2). 

The remaining stars, HD224959, CS22956-028, 
CS22942-019 and LP 625-44, all lie at the short-period 
end of the range produced by our model set A, as is to be 
expected from the above-mentioned selection effects. 

4.2. Model sets G and H: best comparison to observations 

The discussion of the previous section suggests that to bet- 
ter match our models with the observations we should con- 
sider an increase in the amount of third dredge up in low 
mass stars and the effect of switching off thermohaline mix- 
ing. To this end we consider model sets G and H. 
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4.2.1. Model set G: extra third dredge up 

Model sets E, F and G are the same as our default set (set 
A) except that third dredge up is increased in efficiency 
in low mass stars. In model E we have set AM c ^ m ; n = 
— 0.07M Q and A m ; n = 0.8, which is the set of values 
found by Izzard et al. (2004) to be required to match the 
carbon-star luminosity functions of the Magellanic Clouds. 
This results in only a modest increase of the number of 
CEMP stars, see Table 3. A larger effect is obtained by 
setting AM enVjm i n = 0.0 in model F, which allows for 
efficient dredge-up in AGB stars of much smaller initial 
masses and increases the CEMP/EMP ratio to 6.5 per 
cent. Model G is a combination of these parameter choices 
with AM C;m ; n = — 0.1 M Q . With this parameter combina- 
tion all primary TP AGB stars down to initial masses of 
0.8 M Q undergo efficient third dredge up. The IMF peaks 
at low mass, so the number of stars affected is large and 
the CEMP/EMP ratio increases to almost 10% - a fac- 
tor of four increase compared to our default models. The 
NEMP/EMP ratio remains small at 0.3%. 

The distribution of [C/Fe] reaches up to +2.7 and is in 
good agreement with the observations (see Fig. 13a). The 
distributions of nitrogen (Fig. 13a) and sodium suffer the 
same problems as in the default model set. Although the 
total amount of C + N (and hence the total amount of third 
dredge up) as well as the observed trend of [N/Fe] vs [C/Fe] 
roughly match the observations, there is clearly a need for 
some additional CN cycling to convert carbon into nitrogen. 

Model set G has a reduced 13 C efficiency, £13 = 0.1 as 
in model Al. This yields a reasonably good fit to the ob- 
servations for the light-s elements strontium, yttrium and 
zirconium, but the stars most strongly enriched in heavy- 
s elements such as barium (with [Ba/Fe] > +2.7) are not 
well reproduced. Fig. 13b reveals that, firstly, our models 
predict a strong correlation between barium and carbon 
(as well as between other s-process elements) which is not 
seen. The spread in the observations is only partially ex- 
plained by measurement errors. Secondly, the models trace 
out roughly the lowest observed barium abundances for any 
value of [C/Fe]. More barium can be produced in our mod- 
els by further reducing the 13 C efficiency, as discussed in 
Section 4.1.5, but only at the expense of lead. With the 
adopted value of £13 the average observed barium to lead ra- 
tio is well matched, although the most lead-rich stars (with 
[Pb/Fe] > +2.7) are not reproduced. The CEMP-s/CEMP 
ratio for this model set is 94%, in agreement with observa- 
tions. 

Interestingly, some CEMP binaries are made in model 
set G which have periods of one to a few years, similar to 
the observed short-period giant CEMP stars (see Fig. 12). 
These arise from binaries with low-mass AGB primaries 
which have smaller maximum radii and can thus avoid fill- 
ing their Roche lobes in tighter orbits. Because the initial 
mass ratio is close to unity and the primary may even be 
less massive than the secondary due to mass loss prior to 
the AGB, these binaries may also undergo stable RLOF 
without a common envelope phase. 

4.2.2. Model set H: extra third dredge up, no thermohaline 
mixing, common envelope accretion 

Model set H is similar to G, as described in the previous 
section, but is tuned to maximise the CEMP/EMP ratio. 



Thermohaline mixing is turned off so that accreted mate- 
rial remains in the surface convection zone. Before the star 
ascends the giant branch and its convection zone deepens, 
its surface composition is therefore essentially the same as 
that of the primary TPAGB star ejecta. Furthermore, dur- 
ing the common envelope phase accretion of 0.05 M Q of 
material is allowed, so some stars that undergo RLOF can 
become CEMP stars. These changes mean that CEMP stars 
can form out to longer periods and more turn-off stars and 
subgiants, with logg > 3.5, become CEMP stars because 
there is no dilution until first dredge up. While the individ- 
ual effects of these changes on the number of CEMP stars 
are modest (see models C and D in Table 3), in combina- 
tion with extra third dredge-up they yield a CEMP/EMP 
ratio of almost 16%. 

The effect on our model abundance distributions is sig- 
nificant: the sub-giant and turn-off CEMP stars form a 
group of stars with undiluted abundances. In the case of 
[C/Fe] these are between +3 and +4dex (see Fig. 14a), 
which does not match the observations - only one star in 
our observed selection has [C/Fc] = +3 and none are more 
enhanced than this. It may be argued that this excess car- 
bon is converted to nitrogen by extra mixing processes, but 
then [N/Fe] should rise to approximately +3dex in some 
stars, which is not observed either. Our model [N/Fe] distri- 
bution ranges up to +2dex, as do most of the observations, 
but while in our model these correspond to the most C- 
enriched, undiluted stars, the observed nitrogen-rich stars 
have more modest carbon enhancements. 

The undiluted turn-off stars in model H result in a 
peak in the sodium distribution at around 0.8 dex, in better 
agreement with the (broad) peak in the observations than 
our models with thermohaline mixing. Similarly, a broad 
distribution of oxygen abundances is obtained, out to about 
+1.7 dex, also in better agreement with the observations. 

The heavy element distributions are reasonably well 
matched by this model, although the maximum abundances 
predicted (for undiluted stars) are sometimes in excess of 
what is observed. This is especially true for the light-s ele- 
ments, for which this model set makes too many stars with 
large enhancements (up to +2.5 dex for [Zr/Fe] and+2.2 dex 
for [Y, Sr/Fe]) which are not seen in the observations. 
Lanthanum and barium are enriched up to +3.2 dex in this 
model which is in agreement with observations, although 
the most barium-rich objects have too much carbon (see 
Fig. 14b). The models cover the full range of lead observa- 
tions, up to ~ 3.5 dex. 

We must keep in mind that the undiluted stars in our 
models, which have the largest overabundances for all the 
elements discussed above, are also high-gravity stars against 
which there is an apparent observational bias (see Sections 3 
and 4.1.2). They may therefore be overrepresented in our 
model distributions compared to the observational samples. 
We also note that in our observational selection, no corre- 
lation of the abundance distributions against log g is ap- 
parent for any of the elements discussed above. However, 
Denissenkov & Pinsonneault (2008) and Aoki et al. (2008) 
do find evidence for significantly higher average [C/Fe] val- 
ues in turn-off CEMP stars compared to bright giants. 

This model set also includes accretion of up to 0.05 M© 
of material during the common envelope phase. This allows 
CEMP formation in the narrow range of initial separation 
corresponding to primary TPAGB stars which undergo a 
few pulses before overflowing their Roche lobe, and leads 
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Fig. 14. As Fig. 13 for model set H. 



to a peak in the orbital period distribution around one 
year and a tail extending down to 100 days (see Fig. 12). 
This is just the range which includes the stars described in 
Section 4.1.6 and may explain their origin. 

4.3. Model set conclusions 

It appears from the previous sections that our models strug- 
gle to reproduce both the observed frequency of CEMP 
stars and the full range of their abundance patterns. Unless 
we choose a rather extreme combination of model param- 
eters our model CEMP /EMP ratio falls short of the range 
of values deduced from observations and does not exceed 
16% even in our most favourable model. 

In order to come close to reproducing the observations, 
we must assume that TPAGB stars with masses as low 
as 0.8 M experience efficient third dredge up. This yields 
a good fit of the distribution of carbon enhancements, but 
still falls short of reproducing the largest observed s-process 



abundances. If we switch off thermohaline mixing we find, 
on the other hand, that the s-process elements are quite 
well reproduced, but our model CEMP stars have too much 
carbon (and perhaps also lead). We note that both the s- 
process abundance ratios and the CEMP-s/CEMP ratio 
depend on the 13 C pocket efficiency. We can find a reason- 
able match to these constraints by choosing this efficiency 
to be about 0.1 times its default value. 

The orbital periods of CEMP binaries in our observa- 
tional sample agree well with the shortest period CEMP 
stars in both these model sets, so it seems that - at least 
for the few systems with measured periods - the binary 
mass-transfer scenario is compatible with the observations. 

5. Discussion 

We have modelled the observed properties - log g, chemistry 
and orbital period - of CEMP stars and have attempted 
to match our models to observed stars. Our models are 
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successful in reproducing several key observed properties of 
CEMP stars, but struggle to match the full range of the 
observations, most notably the high CEMP fraction, even 
with rather extreme choices of model parameters. 

The first criticism that can be levelled at our work is 
probably our choice of stars from, and subsequent process- 
ing of, the SAGA observational database. We chose a group 
of stars limited by metallicity, [Fc/H] = 2.3 ± 0.5, because 
our stellar models have Z = 1CT 4 . As shown in Section 3, 
the statistics (number of stars, CEMP/EMP ratios etc.) of 
our selection vary remarkably little as a function of [Fe/H] 
so the selection is justified in this respect. However, there 
may be a selection effect inherent in the SAGA database 
because papers are selected for inclusion in it if they contain 
stars with [Fe/H] < —2.5. Stars with higher metallicities, of 
which there are many in the database, are included because 
they just happen to be in those papers. 

We should then consider the selection criterion logg < 
4. While we select identically in the models it is reasonable 
to ask which model constraints are weakened by this choice. 
The answer is few: most stars in the SAGA database are 
giants. This is also the result of a bias against faint, high- 
gravity stars that is inherent in the SAGA database. We 
have made no attempt to model such a selection bias in 
any detail and our simple gravity criterion is not entirely 
successful in selecting against turn-off stars which may be 
affected by the uncertain strength of thermohaline mixing. 
On the other hand, giants are well mixed as their surface 
convection zones deepen on the giant branch and so whether 
we assume efficient thermohaline mixing, or do not, is of 
minor importance provided we look only at giants and avoid 
elements which may be processed in the envelope (see e.g. 
Stancliffe & Glebbeek, 2008 and Stancliffe, 2009). 

The exception to this rule lies with carbon and nitrogen. 
The amount of nitrogen dredged up as the convection zone 
deepens depends on whether accreted carbon has mixed 
deep enough to be burned by the CN cycle. In the case 
of efficient thermohaline mixing this is certainly the case, 
although to some extent the dilution associated with such 
deep mixing reduces the effect of first dredge up. On the 
other hand, if accreted carbon sits at the stellar surface it 
does not burn to nitrogen. Almost all the nitrogen seen in 
the CEMP star must then come from the primary star, pos- 
ing the questions of whether extra mixing, dual-shell flashes 
and dual-core flashes are important. Interesting progress is 
being made regarding these questions (e.g. Stancliffe et al., 
2009). 

We have tried to push the physics of the canonical third 
dredge up as far as possible, inducing it in stars down to 
0.8 M Q with large efficiency. Previous models suggest that 
stars with envelope masses less than 0.5 M do not undergo 
third dredge up, although many of these models are for 
higher metallicity than we are considering here. 

Stancliffe & Glebbeek (2008) have made a model of 
a 0.9 M Q star with Z = 10~ 4 which does undergo third 
dredge up. The efficiency is not high, A = 0.16, but carbon 
is dredged to the surface. At the time that third dredge up 
occurs the star has an envelope mass of just 0.19M Q and 
the dredge up is sufficient to increase the surface carbon 
abundance dramatically, from just 6 x 10~ 6 up to 3.7 x 10~ 3 
by mass fraction. 

If third dredge up is as efficient as assumed in our mod- 
els that provide the best match to observations we would 
expect to see a population of (single) CEMP stars that are 



currently in the TPAGB phase with logg < 0.5. Masseron 
et al. (2006) have suggested that CS 30322-023 is just such 
a star. They construct a 0.8 M©, [Fc/H] = —3.3 detailed 
model sequence with the STAREVOL code (Siess, 2006) for 
comparison with observations of CS 30322-023. They find 
no third dredge up, which is not surprising given that mod- 
els constructed with STAREVOL do not show third dredge up 
unless some kind of neutral-convective-boundary method or 
overshooting is invoked. However the authors make clear 
our lack of quantitative understanding of mixing processes 
in these stars. Other mechanisms such as dual shell and/or 
core flashes, proton ingestion and so-called "canonical" ex- 
tra mixing may, to some extent, mimic the nucleosynthetic 
signature of third dredge up. At the very least, if CS 30322- 
023 is a (single) TPAGB star, as Masseron et al. (2006) sug- 
gest, it certainly has undergone some nucleosynthetic pro- 
cessing in order to reach [C/Fe] ~ +0.6 and [N/Fe] ~ +2.8. 
Convective overshooting may also play a role in enhanc- 
ing third dredge up in low-mass stars: prescriptions such as 
that of Herwig (2000) contain free parameters which have 
the same effect as our AM Cilll i n and A m ; n . 

Even with enhanced third dredge up efficiency our mod- 
els barely make enough CEMP stars to match the ob- 
served fraction. There are more exotic methods for increas- 
ing the number of CEMP stars. One is to accrete mate- 
rial from a TPAGB star on to a main-sequence star during 
the common-envelope phase which should occur for most 
TPAGB primaries that overflow their Roche lobe. We as- 
sumed 0.05 M Q is accreted on to the main-sequence star, 
which is probably too much (Ricker & Taam, 2008), so 
gives an upper limit. In any case, if there is no thermoha- 
line mixing in the main-sequence star even a tiny amount of 
accreted carbon-rich material should turn it into a CEMP 
star. Still, even with 0.05 M Q of accreted material, only a 
few extra CEMP stars are made. This is because the pe- 
riod range in which AGB stars both have enough pulses to 
become carbon rich and then undergo RLOF is rather nar- 
row. At a slightly smaller period RLOF occurs too early, 
i.e. after too few (or no) third dredge up episodes and little 
or no carbon enhancement. 

We also considered a reduction in the common-envelope 
parameter such that acE = 0.1. This does not increase the 
number of CEMP stars - instead it reduces the number of 
non-CEMP stars by forcing many short-period systems to 
merge. After the merger the stellar mass is so large that the 
star evolves quickly and by the present age of the Galaxy it 
is a white dwarf. We do not pretend that this mechanism is 
realistic, but it may help to improve the CEMP/EMP ratio 
match with observations. 

A non-canonical dredge up event may lead to the forma- 
tion of CEMP (and possibly NEMP) stars. Previous works, 
such as Komiya et al. (2007), have suggested that below a 
certain threshold metallicity, in their work [Fc/H] = —2.5, 
dual shell flashes make all the required nitrogen and some 
of the carbon and s-process elements seen in CEMP stars. 
The recent works of Cristallo et al. (2007) and Campbell 
& Lattanzio (2008) indeed show some agreement with the 
observed carbon and nitrogen abundances. However, at the 
metallicity under consideration here ([Fe/H] ~ —2.3) these 
proton-ingestion events are not expected to occur - or 
at least only in the lowest mass stars - and the number 
of CEMP stars should drop as the metallicity increases. 
Extra mixing processes, and a dependence on stellar prop- 
erties other than the metallicity, may well blur the appar- 
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ently sharp boundary between stars that undergo dual-shell 
flashes and those that do not. Incorporation of the results 
of Cristallo et al. (2007) and Campbell & Lattanzio (2008) 
into our models is planned for future work. 

A deliberate omission in the above is discussion of the 
binary fraction in low-metallicity stars. In order to obtain 
anywhere near the observed CEMP fraction we must as- 
sume a 100% binary fraction. This should be compared to 
about 60% in solar-neighourhood G dwarfs (Duquennoy & 
Mayor, 1991) and probably higher among stars more mas- 
sive than the Sun 7 . 

We are left with the situation that in order for our mod- 
els to come close to reproducing the observed CEMP/EMP 
fraction several physical parameters must be pushed to 
the ends of their reasonable range of values. This is not 
a very satisfactory solution, and probably indicates that 
other phenomena require our consideration, such as a shift 
in the initial mass function (Komiya et al., 2007; Lucatello 
et al., 2006), massive-star pollution, primordial supernovae, 
accretion from the interstellar medium etc. All of these so- 
lutions to reproducing the CEMP /EMP fraction have their 
own problems. Most likely, some combination of our physi- 
cal models with, say, a slightly-shifted IMF (or other initial 
distribution), will better reproduce the CEMP/EMP frac- 
tion. Our investigation into this is ongoing and will com- 
prise future work. 

Finally, we note that both the observational statistics 
and the observed abundances of CEMP stars are still uncer- 
tain. Interesting results regarding three-dimensional stellar 
atmosphere models (Collet et al., 2007) may be of relevance 
to CEMP studies. They conclude that the abundances of 
C, N and O in red giants with metallicity [Fe/H] = —3 
may be overestimated by up to 1 dex in traditional one- 
dimensional LTE model atmosphere analyses. If we apply 
these corrections to the Suda database, with the crude as- 
sumption of a linear scaling as a function of metallicity 8 , the 
resulting CEMP/EMP ratio drops to about 14% which well 
fits our enhanced dredge up models. We are not suggest- 
ing this is the answer to the CEMP/EMP ratio problem, 
but it highlights the fact that observed carbon and nitro- 
gen abundances, and hence CEMP number counts, are still 
quite uncertain. 

6. Conclusions 

In an attempt to reproduce the observed CEMP to 
EMP number ratio we have simulated populations of low- 
metallicity ([Fe/H] = —2.3) binary stars with a variety of 
input physics. Our model sets with efficient third dredge up 
in low-mass (down to 0.8 M Q ) stars have CEMP to EMP 
ratios of up to 15%, comparable with the observed ~ 20%. 
They also have low NEMP to EMP number ratios, in agree- 
ment with the observations. Other parameters in our simu- 
lations, such as the efficiency of wind accretion, the common 
envelope parameter etc., have only relatively minor effects 
on our results. 

Acknowledgements. We thank Sara Bisterzo, Simon Campbell, 
Roberto Gallino, Laura Husti, Amanda Karakas, Maria Lugaro and 

7 We completely ignore triples and systems of high multiplicity 
- these are likely to be hierarchical in nature. 

8 The shift in each abundance is then AX dex —fx SX where 
SX is the shift given by Collet et al. (2007) for their [Fe/H] = -3 
models and / = max (0, min [1, {[Fe/H]/ - 3}]). 



the Utrecht stellar evolution group for useful criticism and discussion. 
We thank very much the authors of the SAGA database for their 
willingness to share their database before its publication. RGI thanks 
the NWO for his fellowship in Utrecht and is the recipient of a Marie 
Curie-Intra European Fellowship at ULB. EG acknowledges support 
from the NWO under grant 614.000.303 and NSERC. RJS is funded 
by the Australian Research Council's Discovery Projects scheme un- 
der grant DP0879472. He is grateful to Churchill College for his Junior 
Research Fellowship, under which this work commenced. We would 
also like to thank the referee, Achim Weiss, for many useful sugges- 
tions. 



References 

Abia, C, Recio-Blanco, A., de Laverny, P., et al. 2009, ApJ, 694, 971 
Anders, E. & Grevesse, N. 1989, Geochimica et Cosmochimica Acta, 
53, 197 

Aoki, W., Beers, T. C, Christlieb, N., et al. 2007, ApJ, 655, 492 
Aoki, W., Beers, T. C, Sivarani, T., et al. 2008, ApJ, 678, 1351 
Aoki, W., Norris, J. E., Ryan, S. G., Beers, T. C, & Ando, H. 2000, 
ApJ, 536, L97 

Aoki, W., Norris, J. E., Ryan, S. G., Beers, T. C, & Ando, H. 2002, 
ApJ, 567, 1166 

Barbuy, B., Spite, M., Spite, F., et al. 2005, A&A, 429, 1031 
Beers, T. C. & Christlieb, N. 2005, ARA&A, 43, 531 
Beers, T. C, Preston, G. W., & Shectman, S. A. 1992, AJ, 103, 1987 
Bisterzo, S., Gallino, R., Straniero, O., et al. 2008, in American 

Institute of Physics Conference Series, Vol. 990, First Stars III, 

ed. B. W. O'Shea & A. Heger, 330-332 
BonaCic Marinovic, A., Izzard, R. G., Lugaro, M., & Pols, O. R. 2007, 

A&A, 469, 1013 

Bonacic Marinovic, A. A., Glebbeek, E., & Pols, O. R. 2008, A&A, 
480, 797 

Bondi, H. & Hoyle, F. 1944, MNRAS, 104, 273 

Boothroyd, A. I., Sackmann, I. J., & Ahern, S. C. 1993, ApJ, 416, 762 
Busso, M., Wasserburg, G. J., Nollett, K. M., & Calandra, A. 2007, 
ApJ, 671, 802 

Campbell, S. W. & Lattanzio, J. C. 2008, A&A, 490, 769 
Charbonnel, C. & Zahn, J.-P. 2007, A&A, 467, L15 
Christlieb, N., Green, P. J., Wisotzki, L., & Reimers, D. 2001, A&A, 
375, 366 

Cohen, J. G., Shectman, S., Thompson, I., et al. 2005, ApJ, 633, L109 
Collet, R., Asplund, M., & Trampedach, R. 2007, A&A, 469, 687 
Cristallo, S., Straniero, O., Lederer, M. T., & Aringer, B. 2007, ApJ, 
667, 489 

Denissenkov, P. A. & Pinsonneault, M. 2008, ApJ, 679, 1541 
Dewi, J. D. M. & Tauris, T. M. 2000, A&A, 360, 1043 
Duquennoy, A. & Mayor, M. 1991, A&A, 248, 485 
Eggleton, P. P. 1971, MNRAS, 151, 351 

Eggleton, P. P., Dearborn, D. S. P., & Lattanzio, J. C. 2008, ApJ, 
677, 581 

Eggleton, P. P. & Kiseleva-Eggleton, L. 2002, ApJ, 575, 461 
Frankowski, A. & Jorissen, A. 2007, Baltic Astronomy, 16, 104 
Frebel, A., Christlieb, N., Norris, J. E., et al. 2006, ApJ, 652, 1585 
Fujimoto, M. Y., Iben, I. J., & Hollowell, D. 1990, ApJ, 349, 580 
Gallino, R., Arlandini, C, Busso, M., et al. 1998, ApJ, 497, 388 
Herwig, F. 2000, A&A, 360, 952 
Herwig, F. 2004, ApJ, 605, 425 

Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 
Iben, I. & Renzini, A. 1983, ARA&A, 21, 271 
Iben, Jr., I. 1975, ApJ, 196, 525 

Izzard, R. G., Dray, L. M., Karakas, A. I., Lugaro, M., & Tout, C. A. 

2006, A&A, 460, 565 

Izzard, R. G., Lugaro, M., Karakas, A. I., Iliadis, C, & van Raai, M. 

2007, A&A, 466, 641 

Izzard, R. G. & Tout, C. A. 2004, MNRAS, 350, LI 
Izzard, R. G., Tout, C. A., Karakas, A. I., & Pols, O. R. 2004, MNRAS, 
350, 407 

Johnson, J. A., Herwig, F., Beers, T. C, & Christlieb, N. 2007, ApJ, 
658, 1203 

Jonsell, K., Barklem, P. S., Gustafsson, B., et al. 2006, A&A, 451, 651 
Karakas, A. & Lattanzio, J. C. 2007, Publications of the Astronomical 

Society of Australia, 24, 103 
Karakas, A. I. 2009, MNRAS submitted 

Karakas, A. I., Lattanzio, J. C, & Pols, O. R. 2002, PASA, 19, 515 
Kippenhahn, R., Ruschenplatt, G., & Thomas, H.-C. 1980, 
Astronomy and Astrophysics, 91, 175 



R.G. Izzard et al.: CEMP Binary Populations 



17 



Komiya, Y., Suda, T., Minaguchi, H., et al. 2007, ApJ, 658, 367 
Kroupa, P., Tout, C, & Gilmore, G. 1993, MNRAS, 262, 545 
Lau, H. H. B., Stancliffe, R. J., & Tout, C. A. 2009, MNRAS, 396, 
1046 

Lucatello, S., Beers, T. C., Christlieb, N., et al. 2006, ApJ, 652, L37 
Lucatello, S., Gratton, R. G., Beers, T. C., & Carretta, E. 2005a, ApJ, 
625, 833 

Lucatello, S., Tsangarides, S., Beers, T. C., et al. 2005b, ApJ, 625, 
825 

Luck, R. E. & Bond, H. E. 1991, ApJS, 77, 515 

Lugaro, M., de Mink, S. E., Izzard, R. G., et al. 2008, A&A, 484, L27 
Masseron, T., Johnson, J. A., Plez, B., et al. 2009, A&A accepted, 
ArXiv 0901.4737 

Masseron, T., van Eck, S., Famaey, B., et al. 2006, A&A, 455, 1059 
McClure, R. D. 1984, ApJ, 280, L31 
McClure, R. D. 1997, PASP, 109, 536 

McClure, R. D. & Woodsworth, A. W. 1990, ApJ, 352, 709 
Nelemans, G. & Tout, C. A. 2005, MNRAS, 356, 753 
Nollett, K. M., Busso, M., & Wasserburg, G. J. 2003, ApJ, 582, 1036 
Otsuka, M., Izumiura, H., Tajitsu, A., & Hyung, S. 2008, ApJ, 682, 
L105 

Pols, O. R., Karakas, A. I., Lattanzio, J. C, & Tout, C. A. 2003, 
in Astronomical Society of the Pacific Conference Series, Vol. 303, 
Astronomical Society of the Pacific Conference Series, ed. R. L. M. 
Corradi, J. Mikolajewska, & T. J. Mahoney, 290 

Pols, O. R., Tout, C. A., Eggleton, P. P., & Han, Z. 1995, MNRAS, 
274, 964 

Preston, G. W. & Sneden, C. 2001, AJ, 122, 1545 

Reimers, D. 1975, Circumstellar envelopes and mass loss of red giant 

stars (Springer- Verlag New York, Inc., 1975), 229-256 
Richard, O., Michaud, G., & Richer, J. 2002a, ApJ, 580, 1100 
Richard, O., Michaud, G., Richer, J., et al. 2002b, ApJ, 568, 979 
Ricker, P. M. & Taam, R. E. 2008, ApJ, 672, L41 
Roederer, I. U., Frebel, A., Shetrone, M. D., et al. 2008, ApJ, 679, 

1549 

Ryan, S. G., Aoki, W., Norris, J. E., & Beers, T. C. 2005, ApJ, 635, 
349 

Schlattl, H., Salaris, M., Cassisi, S., & Weiss, A. 2002, A&A, 395, 77 
Schroder, K.-P., Pols, O. R., & Eggleton, P. P. 1997, MNRAS, 285, 
696 

Schuler, S. C, Cunha, K., Smith, V. V., et al. 2007, ApJ, 667, L81 
Siess, L. 2006, A&A, 448, 717 

Sneden, C, Preston, G. W., & Cowan, J. J. 2003, ApJ, 592, 504 
Stancliffe, R. J. 2009, MNRAS, 394, 1051 

Stancliffe, R. J., Church, R. P., Angelou, G. C, & Lattanzio, J. C. 

2009, MNRAS, 396, 2313 
Stancliffe, R. J. & Eldridge, J. J. 2009, MNRAS, 396, 1699 
Stancliffe, R. J. & Glebbeek, E. 2008, MNRAS, 389, 1828 
Stancliffe, R. J., Glebbeek, E., Izzard, R. G., & Pols, O. R. 2007, 

A&A, 464, L57 
Stancliffe, R. J. & Jeffery, C. S. 2007, MNRAS, 375, 1280 
Stancliffe, R. J., Lugaro, M. A., Ugalde, C, et al. 2005, MNRAS, 360, 

375 

Straniero, O., Chieffi, A., Limongi, M., et al. 1997, ApJ, 478, 332 
Suda, T., Katsuta, Y., Yamada, S., et al. 2008, PASJ, 60, 1159 
Thompson, I. B., Ivans, I. I., Bisterzo, S., et al. 2008, ApJ, 677, 556 
Tomkin, J., Lambert, D. L., Edvardsson, B., Gustafsson, B., & Nissen, 

P. E. 1989, A&A, 219, L15 
Tout, C. A. & Eggleton, P. P. 1988, MNRAS, 231, 823 
Tsangarides, S., Ryan, S. G., & Beers, T. C. 2004, MmSAI, 75, 772 
van Loon, J. T., Cioni, M.-R. L., Zijlstra, A. A., & Loup, C. 2005, 

A&A, 438, 273 
Vassiliadis, E. & Wood, P. R. 1993, ApJ, 413, 641 
Weiss, A., Denissenkov, P. A., & Charbonnel, C. 2000, A&A, 356, 181 
Weiss, A. & Ferguson, J. W. 2009, A&A accepted, ArXiv 0903.2155 
Weiss, A., Schlattl, H., Salaris, M., & Cassisi, S. 2004, A&A, 422, 217 



Appendix A: Dredge Up Prescriptions 

A.l. First dredge up 

The change in surface abundance of isotopes j at first 
dredge up, AX,, is interpolated from a table of detailed 
models with Z = 1CP 4 in mass range 0.5 < M/M Q < 12. 
A correction factor / C no = X C no(TMS)/X C no(ZAMS), 



the ratio of CNO mass fraction at the terminal-age main 
sequence (TMS) and zero-age main sequence (ZAMS), is 
then applied to CNO elements to take into account accre- 
tion during the main sequence. 

In the Izzard et al. (2006) model first dredge up is 
considered as an instantaneous event. In terms of time 
evolution this is a reasonable assumption because giant- 
branch evolution is fast, but in terms of luminosity or 
gravity this approximation is not good and it proves 
difficult to compare to e.g. the [C/Fe] vs log(L/L ) 
data of Lucatello et al. (2006). To resolve this problem 
the changes in abundances are modulated by a factor 
/ p = min [(M c - M CiB agb) / (^c.idupmax - M^bagb) , 1] 
where M c is the core mass, Mc,idupmax is the core mass 
at which first dredge up reaches its maximum depth and 
-Mc.bagb is the core mass at the base of the giant branch, 
before first dredge up starts. M Cj bagb is known from the 
stellar evolution prescription and Mc,idupmax, is interpo- 
lated from a grid of models constructed with the TWIN stel- 
lar evolution code (Eggleton & Kiseleva-Eggleton, 2002). 

In summary, the surface abundances changes at first 
dredge up are given by /cNo/pAXj. They agree well with 
the detailed models, as a function of M c , logL and time. 

A. 2. Third dredge up 

Abundance changes at third dredge up are treated in a 
similar way to the prescription of Izzard et al. (2004) and 
Izzard et al. (2006). Intershell abundances are interpolated 
from tables based on the Karakas et al. (2002) detailed 
models the metallicities of which extend down to Z = 10 . 

In low-metallicity TPAGB stars dredge up of the 
hydrogen-burning shell enhances the surface abundance of 
13 C and 14 N (at higher metallicity the effect is negligible 
because the initial abundance of 13 C and 14 N is relatively 
large). This is modelled by dredging up SM of hydrogen- 
burnt material during each third dredge up, where the 
abundance mixture in this material is enhanced in 13 C and 
14 N according to 



X C i3 = 0.006 x A C i2 and 
Ani4 = 0.28 x Xci2 , 



(A.l) 
(A.2) 




x min 1, 



A T P 



1 2 N 



10 



(A.3) 



and M(t) is the instantaneous stellar mass, M cnv (t) is the 
instantaneous envelope mass, Ntp is the thermal pulse 
number, X12 is the envelope abundance of 12 C and e = 
2Q-20 rj ne £ rg ^. ^ erm gj ves the amount of H-burnt mate- 
rial dredged up, the second term is a turn-on effect as the 
star reaches the asymptotic regime and the third term is a 
turn-off effect for small envelopes. 



Appendix B: Mass-loss prescriptions 

We consider three mass-loss prescriptions for TPAGB stars. 

VW93 The formalism of Vassiliadis & Wood (1993, VW93) 
relates the mass-loss rate to the Mira pulsation period 
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of the star, given by 

p Q = 1O -2.O7-O.9OL/L +1.941og lo (ii/R ) days ( R 

The mass loss rate is then given by, as in Karakas et al., 
2002, i.e. without the M/M Q - 2.5 term of the original 
VW93 prescription, 

M = Mywo = /vw (-H.4 + 0.0125P ) M Q yr" 1 (B.2) 

unless Pq > (500 — APvw) days in which case a super- 
wind is applied 



M = max (Mvwo./vw— ) 



(B.3) 



where 

v w = 10 5 (-13.5 + 0.056P ) cms^ 1 . (B.4) 

The free parameters /vw and APyw subtly affect the 
mass-loss rate. The factor /vw is a simple multiplier, 
which is 1 by default (see model set 27). The period 
shift APvw allows the onset of the superwind to be 
delayed, e.g. APyw = —100 days in model set 33 - it is 
zero by default. 
Reimers The Reimers mass-loss rate is given by 



M = 4 x 10" 13 r/ 



j| ^M 9 yr- 

M UJ 



(B.5) 



where 77 is a parameter of order unity (Reimers, 1975) 
which we vary in model sets 10, 11 and 12. 
van Loon In model set 13 we use the split form of van Loon 
et al. (2005) appropriate to oxygen-rich red giants, 



log 10 [MAMoyr- 1 ) 



{ 



-5.6+1.10Z 4 
-5.3 + 0.82/4- 



- 5.2t 35 l 4 < 0.9 and 
10.8% Z 4 >0.9, 



(B.6) 



where I4 = log 10 (i/10 4 L Q ) and £35 = 
log 10 (T cff /3500K). Note, if T cS > 4000 K we en- 
force a minimum mass-loss rate of 10~ 4 M Q year -1 
because the above formula can approach zero as the 
temperature rises (and the envelope mass becomes 
small) as a star approaches the white-dwarf cooling 
track. 



Appendix D: All Model Sets and Results 

Table D.l shows the full set of models we considered, of 
which Table 1 is a subset. 

Table D.2 shows the full set of CEMP, CNEMP and 
NEMP to EMP ratios for all our model sets of which Table 3 
is a subset. 



Appendix E: Observation database 

Our observational selection is taken from the SAGA 
database as compiled by Suda et al. (2008) combined with 
data Lucatello et al. (2006). 

When data exists for the same star from more than one 
source, we take the arithmetic mean of the values and add 
errors in quadrature. In the case of log-values, e.g. [Fe/H] or 
logg, we simply average the log- values rather than attempt 
a more sophisticated approach. This makes little difference 
to our final results. In the case of data limits (e.g. x < 4) we 
ignore the data - few data are of this type and the general 
result is not affected. 

We ignore error bars in the sense that, e.g. a star with 
[Fe/H] = —2.9 ± 0.2 is not included in our selection, even 
though it may well have - in reality - [Fe/H] = — 2.7 and 
hence qualify. This is the price we pay for a simple selection 
procedure and in the large number limit (the database has 
about 1300 stars) it is not a problem. 



Appendix C: Binary distributions 

Our default binary-star distribution is the combination of 

1. The initial mass function (IMF) of Kroupa et al. (1993, 
KTG93) for the initial primary mass Mi 

Mi/M < m 



V(Mi) = ^ 



ai(Mi/M Q ) P1 m < M 1 /M Q < m 



a 2 (Mi/M Q ) P2 mi < Mi/M < m 2 
a 3 (Mi/M Q ) P3 m 2 < Mi/M Q < 
m > m max 

(C.l) 

where p\ — —1.3, p 2 = —2.2, ps = —2.7, m,Q = 0.1, 
mi = 0.5, m-2 = 1.0 and m max = 80.0. Continuity re- 
quirements and J ip(M)dM = 1 give the constants en, 
<Z2 and 03. 

2. A distribution flat in q = M2/M1 for the initial sec- 
ondary mass M2, where M 2 < Mi 

3. A distribution flat in In a (i.e. probability ~ 1/a) for the 
separation a where 3 < a < 10 5 . 

4. Initially circular binaries (except for model set Ae5). 
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lVJ.OU.ei Oct 


Physical Parameters (differences from model set A^j 


A 

A 

/~i 771 -f rtJ?Q 


q:ce =0.1 and respectively 


D 

JD 


Bondi-Hoyle cubh — 5 


Situ 

Q 
O 


Initial eccentricity 6 — 0.5 
v^ii/\r^ parameter jd — iu 


E 1 


Imra JJUF calibration: i\M c , m i n = — U.UY Mq, A m i n = U.o 


1(1 11 10 

1U, 1 1 , 1 z 


Reimers AGB wind: ?y = 0.1, 1 and 5 respectively 


1 Q 


AGB wind of Van Loon 


11 1 K H 
14, -*<->, O 


Common envelope accretion 0.01, 0.02 and 0.05 Mq respectively 


Al, A2 


C efficiency parameter 0.1 and 0.01 respectively 




Thermohaline mixing disabled 


A 4- 1 <n A-t-Q 
AtlZ, Ato 


Minimum CEMP age 12 and 8 Gyr respectively 


TTl 
_^ 


Imrd DUr -Menv.min = OMq 


Z 1 


M e nv,min = UMq, $13 = 0.01, AiW Cjmin = —U.l Mq, Amin = (J. 5, /vw = U.l 


28 


71/T n 1\/T /- A Al A 71 /T A 1 AT \ A K 
Menv.min = M Q , $13 = 0.01, /\M c ,mm = — 0.1 Mq, A m in = 0.5 


Apo, Ap / 


[C/Fe] min — 0.5 and 0.7 respectively 


31, 32 


As 27 with [C/Fe] min = 0.5 and 0.7 respectively 


33 


A n i riA j „ 

Arvw = — lUUdays 


34, Bl 


71 /f A T\ T r 1 " A Al A 71 /f A 1 A 4" Al J A O 

Menv.min = Mq , $13 = (J. 01, AAi c . m in = —0.1 M© with A m i n = 0.1 and 0.8 respectively 


36 


Menv.min = 0.25 M© , $13 = 0.01, AMc.min = — 0.1 Mq, A m i n = 0.8 


37, 38, 39, 4O 


"NT 1 1 ' J.' ACT11CTJO 

JNelemans common-envelope prescription, 7 = 0.5, 1, 1.5 and 2 respectively 


41 

42 


As but Nelemans prescription lor <? > 0.2 only (otherwise the delault prescription) 
As but Nelemans prescription only for the first common-envelope phase 


D ft 

BZ 


Qce = 0.1, Menv.min = UMq, AiMt c . m in = —U.l Mq , A m in = O.o , uomenv accretion 0.00, INo tnermonalme 


44 


38 and B2 combined 


45 


As Bl with $13 = 0.1 


48 


71/T AA/T A71/T A "1 A If \ ACr 1 " AA1 
Menv.min = M© , AM c , m i n = —0.1 M© , A min = 0.5, $13 = 0.01 


G 


Menv.min = M© , AM c , m in = -0.1 Mq, A min = 0.8, $13 = 0.1 


50 


Menv.min = Mq , AM c . m i n = -0.1 Mq, Amin = 0.5, $13 = 0.001 


51, 52, 53 


Comenv accretion 0.05 Mq, Menv.min = OMq, 

AMe.min = -0.1 Mq, Amin = 0.5, 

$13 = 0.1, 0.01 and 0.001 respectively 




Pntnenv arrrpfinn 0^ A/T^ fl/f — ft l\ I 
V^UIIlcIlv cLIUll U.UtJIVlQ, JvJenv,min — U 1V10 , 

AMe.min = -0.1 Mq, Amin = 0.8, 

$13 = 0.1, 0.01 and 0.001 respectively 


57, 58, 59 


Menv.min = Mq , AM c , m i n = -0.1 Mq, 
Amin = 0.5, 

$13 = 0.1, 0.01 and 0.001 respectively 



Table D.l. The full list of our binary population models (a subset is shown in Table 1). The meanings of the symbols are given 
in Section 2. 
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Model Set 


CEMP/EMP 70 


CJNhjMP/lijMP 70 


"VTT - 'A T T> i'T~ ''TV T T> 0~/ 

JNBjMP/BjMP 7o 


A, Al, A2 


inn 1 a ao 1 

2.300 ± 0.034 


A AAO 1 A AAA 

0.098 ± 0.002 


A .~>/"'^7 1 A A A"" 

0.267 ± 0.007 


CEpl 


a t a 1 a nnA 

2.370 ± 0.029 


A 1 AA 1 A AAA 

0.100 ± 0.002 


A A TO 1 A AAT 

0.273 ± 0.007 


CE3 


A OAA 1 A AOO 

2.300 ± 0.033 


A AAT 1 A AAA 

0.097 ± 0.002 


A A/ 1 T 1 A AA/ 1 

0.267 ± 0.006 


B 


nrn 1 a a 

3.050 ± 0.038 


A 1 f IT 1 A AAO 

0.157 ± 0.003 


A OAO 1 A AA/ 1 

0.323 ± 0.006 


Ae5 


a -1 r a 1 a aoo 

2.150 ± 0.033 


A AO/ 1 1 A AAA 

0.086 ± 0.002 


A A 1 A I A AAT 

0.249 ± 0.007 


8 


O A cn 1 A AO /I 

2.460 ± 0.034 


A 1 A 1 1 A AAA 

0.101 ± 0.002 


A A£?0 I A C\C\G 

0.268 ± 0.006 


E 


a aaa 1 a n <n 

2.900 ± 0.040 


A AAO 1 A AAA 

0.098 ± 0.002 


A A/" 1 T 1 A AA/ 1 

0.267 ± 0.006 


10 


A /"* /) A 1 a nnn 

2.640 ± 0.032 


A AAO 1 A AAO 

0.208 ± 0.003 


A O 1 A 1 A AAI - 

0.310 ± 0.005 


11 


1 Ar\f\ 1 a a a -1 

1.490 ± 0.021 


A A f O 1 A AA1 

0.053 ± 0.001 


A AAI 1 A AAA 

0.091 ± 0.002 


12 


/-i a oa 1 a a a a 

0.489 ± 0.009 


A AAI"" 1 A AAA 

0.005 ± 0.000 


A Al t~ 1 A AAI 

0.015 ± 0.001 


13 


a 1 ai 1 a aa a 

0.101 ± 0.004 


0.000 


0.000 


U 


a ,1 /"* a 1 a nor 

2.460 ± 0.035 


A AAA 1 A AAA 

0.099 ± 0.002 


A AAI 1 A AAT 

0.291 ± 0.007 


15 


a /> <a 1 a ao/ 1 

2.640 ± 0.036 


A 1 AA 1 A AAA 

0.100 ± 0.002 


A OAO 1 A AAT 

0.303 ± 0.007 


C 


a a a 1 a aoo 

2.940 ± 0.038 


A 1 AO 1 A AAA 

0.103 ± 0.002 


A 1 1 1 A AAT 

0.311 ± 0.007 


D 


,i 01 a 1 a a a 

4.210 ± 0.039 


A AAA 1 A AAO 

0.290 ± 0.003 


A A AA 1 A AA A 

0.409 ± 0.004 


Atl2 


aoa 1 a a 

2.230 ± 0.047 


A AAO 1 A AAO 

0.093 ± 0.003 


A A t~ O 1 A AAA 

0.258 ± 0.009 


At8 


a /i oa 1 a a a 

2.480 ± 0.028 


A 1 1 A 1 A AAA 

0.110 ± 0.002 


A ATO 1 A AAI - 

0.278 ± 0.005 


F 


6.470 ± 0.030 


A 1 AO 1 A AAA 

0.103 ± 0.002 


A A J ' *~ 1 A AA/ 1 

0.267 ± 0.006 


21 


O O ,1 A 1 A A /I 1 

8.840 ± 0.041 


A 1 1 A AAO 

0.171 ± 0.003 


A ATO 1 A AA A 

0.273 ± 0.004 


28 


O OAA 1 A AOA 

8.320 ± 0.039 


A 1 AO 1 A AAA 

0.103 ± 0.002 


A A/ 1 / 1 1 A AA/ 1 

0.266 ± 0.006 


A r 

Apo 


O A £?A 1 A A /I O 


A AAO 1 A AAO 

0.203 ± 0.003 


A A^T 1 A AA A 

0.267 ± 0.004 


A rt 

Ap7 


O AAA 1 A AOA 

2.990 ± 0.039 


A 1 CT O _l_ A AAO 

0.158 ± 0.003 


A A^T 1 A AACT 

0.267 ± 0.005 


31 


1 1 CAA 1 A ACQ 

11.600 ± 0.053 


A Attr* 1 A AAO 

0.256 ± 0.003 


A ATO 1 A AA A 

0.273 ± 0.004 


32 


1 A A AA 1 A A A O 

10.400 ± 0.048 


A AAI 1 A AAO 

0.221 ± 0.003 


A 070 1 A AA A 

0.273 ± 0.004 


33 


O A AA 1 A AOCT 

2.490 ± 0.035 


A 1 1 G _l_ A AAO 

0.116 ± 0.003 


A O A A 1 A AAO 

0.340 ± 0.008 


34 


/1 r AA 1 A AOA 

6.520 ± 0.030 


A 1 AO 1 A AAA 

0.103 ± 0.002 


A A/ 1 T 1 A AA/ 1 

0.267 ± 0.006 


Bl 


A A OA 1 A A 1 1 

9.430 ± 0.044 


A 1 AO 1 A AAA 

0.103 ± 0.002 


A A/ 1 / 1 1 A AA/ 1 

0.266 ± 0.006 


36 


CT C A A I A A£?A 

5.640 ± 0.060 


A 1 AA 1 A AAA 

0.102 zb 0.002 


A A^T 1 A AA£? 

0.267 ± 0.006 


37 


A A /"* A I A AO/ 1 

2.260 ± 0.036 


A AA/ 1 1 A AAA 

0.096 ± 0.002 


A A/ 1 A 1 A AA/ 1 

0.262 ± 0.006 


38 


A A "A 1 A AO/ 1 

2.270 ± 0.036 


A AA/ 1 1 A AAA 

0.096 ± 0.002 


A A/ 1 A 1 A AA/ 1 

0.264 ± 0.006 


39 


A AOA 1 A A O f 

2.280 ± 0.035 


A AAT 1 A AAA 

0.097 ± 0.002 


A A/ 1 f 1 A AAT 

0.265 ± 0.007 


40 


A OAA 1 A AO/ 1 

2.300 ± 0.036 


A AAO 1 A AAA 

0.098 ± 0.002 


A A/ 1 r" 1 A AAT 

0.265 ± 0.007 


41,42 


A AOA 1 A AOI - 

2.280 ± 0.035 


A AAT 1 A AAA 

0.097 ± 0.002 


A A/ 1 f 1 A AA/ 1 

0.265 ± 0.006 


B2 


1 A AAA 1 A A /"* O 

14.900 ± 0.063 


A OAO 1 A AAO 

0.303 ± 0.003 


A (fir 1 A A A 1 

0.425 ± 0.004 


44 


1 1 S~' A A 1 A A/ 1 / 1 

14.600 ± 0.066 


A AO/ 1 1 A AAO 

0.286 ± 0.003 


A OOA I A AA A 

0.389 ± 0.004 


40, Lr 






n 966 4- n nn6 


47 


2.060 ± 0.028 


0.000 


0.002 ±0.000 


48, 50 


8.320 ± 0.039 


0.103 ±0.002 


0.266 ±0.006 


51, 52, 53 


14.700 ± 0.064 


0.298 ± 0.003 


0.426 ± 0.004 


H, 55, 56 


15.500 ±0.068 


0.298 ± 0.003 


0.426 ± 0.004 


57, 58, 59 


12.900 ±0.057 


0.278 ± 0.003 


0.387 ±0.004 



Table D.2. Percentage of CEMP, CNEMP and NEMP (sub-)giants relative to total EMP giants in all our model binary populations 
(see Section 2.6 for selection criteria). The errors convey Poisson statistics only. The final column gives the number ratio of CEMP-s 
to CEMP stars. 



